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ABSTRACT 

A  thermistor  bolometer  has  been  used  to  determine  the  radiant 
energy  transfer  from  a  typical  laboratory  plasma.  Radiation  has  been 
found  to  be  a  significant  mechanism  of  energy  transfer  from  an  atmos¬ 
pheric  argon  plasma  jet,  amounting  to  as  much  as  18  percent  of  the 
net  power  input  to  the  plasma.  The  radiant  energy  transfer  from  an 
atmospheric  nitrogen  plasma  jet  under  comparable  conditions  has  been 
found  to  be  less  than  1  percent  of  the  net  power  input  to  the  plasma. 
The  radiative  power  density  in  an  atmospheric  argon  plasma  h  £i  s  Id  ©  g  n 
determined.  Contributions  from  both  bremsstrahlung  and  recombination 
radiation  have  been  found  to  be  important.  A  model,  based  on  the 
Kramers  and  the  UnsHld  approximations,  has  been  used  to  predict  the 
radiative  power  density  in  an  atmospheric  argon  plasma. 
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CHAPTER  1 

INTRODUCTION 

1.1  IMPORTANCE  OF  RADIATIVE  PROCESSES  IN  GAS  DYNAMICS 

Work  in  the  field  of  gas  dynamics  has  frequently  ignored  radiative 
processes*  The  justification  for  this  approximation  is  that  most  gas 
dynamics  problems  are  concerned  with  relatively  cool  gases  in  which 
radiative  processes  are  insignificant*  In  the  realm  of  high  temperature 
gas  dynamics,  however,  radiation  must  be  considered  as  a  potentially 
significant  mechanism  of  energy  transport. 

With  the  growing  interest  in  flow  problems  involving  plasmas, 
plasma  generating  facilities  are  becoming  familiar  items  of  laboratory 
equipment.  The  plasma  torch  is  a  common  method  of  generating  plasmas 
in  the  laboratory*  The  subject  of  this  report  is  the  measurement  of 
the  radiative  energy  transfer  from  plasmas  generated  in  a  typical 
laboratory  plasma  torch. 

1.2  CLASSICAL  APPROACH  TO  RADIATIVE  HEAT  TRANSFER 

The  basic  relationship  in  problems  involving  radiative  heat  trans¬ 
fer  is  the  Stefan-Boltzmann  law,  which  relates  the  emissive  power  of  an 
ideally  black  surface  to  its  temperature 

WB  =  <r  T4  (1-D 

2 

where  Wg  is  the  emissive  power  of  an  ideally  black  surface  (Btu/hr  ft  ) 
<r  is  the  Stefan-Boltzmann  constant,  0.1713  x  10 Btu/hr  ft2  °R^ 

T  is  the  absolute  temperature  of  the  surface  (°R) 
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The  spectral  distribution  of  black-body  radiation  is  given  by  Planck's 
law 


2trhc2 

Xs  (tfCh/kXT  _  |) 


(1-2) 


where  is  the  monochromatic  emissive  power  of  a  black  surface  at 

q 

wavelength  X  (ergs/sec  cur  ) 

*»27 

h  is  Planck's  constant,  6.63  x  10  erg  sec 

xo 

c  is  the  velocity  of  light,  3«0  x  10  cm/sec 
X  is  the  wavelength  of  radiation  (cm) 
k  is  Boltzmann's  constant,  1«38  x  10”^  erg/°K 
T  is  the  absolute  temperature  of  the  surface  (°K) 

The  emissive  power  of  a  real  surface  is  always  less  than  the 
emissive  pother  of  a  black  surface  at  the  same  temperature.  The  classical 
approach  to  the  problem  of  radiative  heat  transfer  between  solid  surfaces 
uses  the  definition  of  the  emissivity  of  a  real  surface 


W_ 

**  WB 


(1-3) 


where  e  is  the  emissivity  of  the  real  surface 

2 

W  is  the  emissive  power  of  the  real  surface  (Btu/hr  ft  ) 

Wg  is  the  emissive  power  of  a  black  surface  at  the  same  temperature 
(Btu/hr  ft2) 

In  general,  the  emissivity  of  a  real  surface  is  a  function  of  wavelength 
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and  temperature .  If  the  radiation  emitted  by  a  particular  surface  has 
the  same  spectral  distribution  as  the  radiation  from  a  black  surface  at 
the  same  temperature,  the  emissivity  of  that  surface  is  independent  of 
wavelength.  A  surface  having  this  property  is  referred  to  as  a  gray 
body.  For  most  engineering  purposes,  radiators  are  considered  as  gray 
bodies.  In  this  case,  the  emissive  power  of  a  real  surface  is  expressed 
as 

W  —  e  o-T4  (1-k) 

2 

where  ¥  is  the  emissive  power  of  the  real  surface  (Btu/hr  ft  ) 

€  is  the  emissivity  of  the  surface 

<y  is  the  Stefan-Boltzmann  constant,  0.1713  x  10*"^  Btu/hr  ft^  °R^ 

T  is  the  absolute  temperature  of  the  surface  (°R) 

¥hereas  the  proportionality  between  the  emissive  power  of  a  black  surface 
and  the  fourth  power  of  its  absolute  temperature  can  be  derived  from 
thermodynamic  considerations,  Eq.  ±-h  is  an  empirical  relationship. 

The  spectrum  of  radiation  emitted  from  a  plasma  does  not  resemble 
black-body  radiation.  Since  a  plasma  is  nearly  transparent,  the  radia¬ 
tion  is  emitted  throughout  the  volume  of  the  plasma.  Very  strong  tempera¬ 
ture  gradients  exist  in  a  plasma  jet.  For  these  reasons,  the  classical 
approach  to  radiative  heat  transfer  embodied  in  Eq.  1-1|  is  an  inappro¬ 
priate  and  inconvenient  method  to  describe  the  radiation  emitted  from  a 
plasma  jet. 

1.3  MACROSCOPIC  DESCRIPTION  OF  RADIATION  EMITTED  FROM  A  PLASMA  JET 

The  radiation  emitted  from  a  transparent  plasma  is  a  unique  function 
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of  the  thermodynamic  state  of  the  plasma  provided  the  plasma  is  in 
local  thermodynamic  equilibrium®  In  order  to  calculate  the  radiant 
energy  emitted  from  a  transparent  plasma  of  any  given  geometrical  con¬ 
figuration,,  the  thermodynamic  state  of  the  plasma  at  every  point  and  the 
radiation  emitted  per  unit  volume  of  the  plasma  as  a  function  of  the 
thermodynamic  state  of  the  plasma  must  be  known® 

In  determining  the  radiation  emitted  from  a  plasma  jet  generated 
with  a  plasma  torch*  this  approach  meets  with  a  serious  practical  diffi¬ 
culty®  For  a  plasma  jet  issuing  into  atmospheric  pressure*  a  knowledge 
of  the  thenaodynamic  state  of  the  plasma  is  provided  by  a  knowledge  of 
the  temperature  at  every  point  in  the  plasma  jet®  The  determination  of 
the  temperature  at  eveiy  point  in  the  plasma  jet  is  a  difficult  task. 

It  is  desirable  to  find  some  more  easily  determined  property  of  the 
plasma  jet  to  use  as  an  independent  variable  when  measuring  the  radia¬ 
tion  emitted  from  the  plasma  jet. 

The  net  power  input  to  the  plasma  is  used  in  this  report  as  the 
independent  variable  with  which  the  macroscopic  radiation  measurements 
are  correlated.  The  use  of  the  net  power  input  in  this  context  requires 
the  implicit  assumption  that  the  temperature  distribution  in  the  plasma 
jet  at  a  given  flow  rate  is  a  unique  function  of  the  net  power  input. 

The  validity  of  this  assumption  for  the  plasm  jets  studied  in  this 
investigation  has  been  verified  indirectly. 

1.4  MICROSCOPIC  DESCRIPTION  OF  RADIATION  EMITTED  FROM  A  PLASMA 

As  discussed  in  Section  1®3$  the  radiation  emitted  per  unit  volume 
from  a  plasma  in  local  thermodynamic  equilibrium  is  a  unique  function 
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of  the  thermodynamic  state  of  the  plasma.  This  relationship  has  been 
determined  experimentally  for  an  argon  plasma  at  atmospheric  pressure. 
These  results  are  only  applicable  to  transparent  argon  plasmas  at  at¬ 
mospheric  pressure  which  are  in  local  thermodynamic  equilibrium,  but 
are  not  otherwise  restricted. 

1.5  EXPERIMENTAL  INVESTIGATIONS 

The  intensity  of  radiation  emitted  from  a  plasma  jet  generated  in 
a  plasma  torch  and  exhausted  to  atmospheric  pressure  has  been  measured 
with  a  thermistor  bolometer.  The  results  presented  in  Chapter  3  are 
based  on  the  macroscopic  description  of  the  plasma  jet  using  the  net 
power  input  as  the  independent  variable.  The  results  include  determina¬ 
tions  of  the  total  radiation  from  an  argon  plasma  jet,  the  total  radia¬ 
tion  from  a  nitrogen  plasma  jet,  and  the  total  radiation  from  an  argon 
plasma  jet  seeded  with  carbon  particles.  Pointwise  radiation  measurements 
were  used  to  determine  the  distribution  of  the  radiative  flux  along  the 
axis  of  the  argon  plasma  jet.  The  angular  distribution  of  the  radiation 
from  the  argon  plasma  jet  and  the  plasma  torch  was  measured  and  was  used 
to  calculate  the  actual  radiative  heat  transfer  from  the  plasma  jet  and 
torch.  Measurements  of  the  temperature  and  radiation  distributions 
leading  to  a  determination  of  the  intensity  of  radiation  emitted  per 
unit  volume  from  an  atmospheric  argon  plasma  as  a  function  of  tempera¬ 
ture  are  discussed  in  Chapter  £. 
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CHAPTER  2 

EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

Detailed  descriptions  of  the  plasma  generating  facility  and  the 
experimental  procedure  are  given  in  Appendix  A. 

2.1  PLASMA  TORCH 

The  plasmas  studied  in  these  experiments  were  generated  in  a 
Model  F-1;0  plasma  torch  manufactured  by  Thermal  Dynamics  Corporation. 

In  the  plasma  torch  the  gas  is  heated  as  it  flows  through  a  direct 
current  arc  between  two  electrodes,  thereby  generating  the  plasma.  The 
gas  enters  the  arc  chamber  through  a  straight-flow  ceramic  disc.  The 
replaceable  copper  anode  acts  as  the  exit  nozzle  for  the  plasma.  The 
replaceable  tungsten  cathode  is  located  upstream  from  the  anode.  The 
arc  is  formed  in  the  annular  region  between  the  nozzle  (anode)  and  the 
electrode  (cathode).  A  #57  x  15  nozzle  and  a  #i|-A  electrode  were  used 
to  generate  an  argon  plasma j  a  #57  x  22  nozzle  and  a  #3~M  electrode 
were  used  to  generate  a  nitrogen  plasma.  Both  the  electrode  and  nozzle 
are  water  cooled. 

2.2  BOLOMETER 

A  Model  1312  thermistor  bolometer  manufactured  by  Servo  Corpora¬ 
tion  of  America  was  used  to  detect  the  radiation  from  the  plasma.  The 
bolometer  sensing  element  is  a  small  flake  of  semiconductor  material, 

0.5  nroi  square,  which  has  a  large  negative  temperature  coefficient  of 
resistance.  The  flake  is  blackened  so  that  it  absorbs  about  80  per  cent 
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of  the  radiation  incident  upon  it  (Ref.  l). 

The  bolometer  sensing  element  is  mounted  in  a  sealed  cylindrical 
case  5/8  in.  long  by  5/8  in.  in  diameter  as  shown  in  Figure  2-1. 
Radiation  is  transmitted  to  the  bolometer  sensing  element  through  a 
sodium  chloride  window  which  transmits  radiation  in  the  frequency  range 
from  0*175  x  lO^4  sec'”'1'  to  11*2  x  10"^  sec"'1'  (Ref.  2).  A  second  elements 
shielded  from  radiation  but  otherwise  identical  to  the  bolometer  sensing 
element,  is  also  mounted  in  the  cylindrical  case.  This  element  is  used 
to  compensate  for  variations  in  the  ambient  temperature. 

Both  bolometer  elements  are  part  of  the  bridge  circuit  shown  in 
Figure  2-2.  A  Ballantine  Model  801  a.c.  voltmeter  is  used  to  measure 
the  output  signal  of  the  bridge  amplifier.  The  radiation  incident  on 
the  bolometer  is  chopped  at  30  c.p.s.  to  provide  an  a.c.  signal  to  the 
amplifier.  The  chopping,  is  accomplished  with  a  6  in.  disc  having  30 
equally  spaced  1/4  in.  holes  which  is  rotated  at  1  r.p.s*  by  a  synchro¬ 
nous  motor.  The  chopper  assembly  is  mounted  on  top  of  the  box  contain¬ 
ing  the  bolometer,  bridge,  and  amplifier*  The  power  supplies  are 
external . 

2.3  EXPERIMENTAL  ARRANGEMENT 

A  photograph  of  the  basic  arrangement  of  the  experimental  apparatus 
is  shown  in  Figure  2-3.  The  bolometer  unit  is  mounted  on  a  laboratory 
jack  which  permits  adjustment  of  its  height.  The  laboratory  jack  is 
mounted  on  a  table  which  can  be  moved  in  two  horizontal  directions  with 
lead  screws.  This  table  rests  on  rails  which  run  the  length  of  the 
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Fig.  2-3  Photograph  of  Experimental  Apparatus 
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bench.  With  the  unit  clamped  to  the  rails  at  some  nominal  location, 
the  bolometer  sensing  element  can  be  moved  about  nine  inches  in  any 
direction  in  space.  The  position  of  the  sensing  element  is  determined 
by  indicating  dials  on  the  lead  screws  and  by  a  vertical  scale  attached 
to  the  laboratory  jack. 

The  plasma  torch  is  operated  from  a  control  panel.  The  current 
to  the  plasma  torch  and  the  temperature  rise  of  the  cooling  water  were 
recorded  continuously.  Instruments  indicating  the  temperature,  pressure, 
and  flow  rate  of  the  gas  and  the  flow  rate  of  the  cooling  water  are 
housed  in  the  control  panel.  The  voltmeter  indicating  the  voltage  drop 
across  the  plasma  torch  and  the  voltmeter  indicating  the  bolometer  out¬ 
put  are  located  on  the  bench. 

2.4  EXPERIMENTAL  PROCEDURE 

The  output  signal  of  the  bolometer  is  proportional  to  the  radia¬ 
tive  flux  incident  on  the  bolometer  element,  as  determined  by  calibrating 
the  instrument.  With  appropriate  mathematical  models  and  physical 
arrangements  of  the  apparatus,  the  bolometer  was  used  to  measure  the 
total  radiation  from  the  plasma  jet,  the  actual  radiative  heat  transfer 
from  the  plasma  jet  and  the  plasma  torch,  the  distribution  of  the  radia¬ 
tive  flux  along  the  axis  of  the  plasma  jet,  and  the  internal  radiative 
power  density  in  the  plasma  jet.  Since  the  net  power  input  to  the 
plasma  jet  is  used  as  the  independent  variable  with  which  to  represent 
the  macroscopic  radiation  measurements  and  as  the  parameter  of  torch 
operation  with  the  microscopic  radiation  measurements,  the  plasma  torch 
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is  instrumented  to  provide  knowledge  of  the  instantaneous  net  power 
input  to  the  plasma.  During  the  experimental  runs  the  bolometer  output 
signal  and  the  voltage  drop  across  the  plasma  torch  were  recorded  simul¬ 
taneously,  noting  the  time  of  each  reading  on  a  stopwatch  synchronized 
with  the  current  and  the  cooling  water  temperature  recorders. 

Descriptions  of  the  physical  arrangements  of  the  bolometer  and 
plasma  torch  are  given  in  the  sections  of  the  report  describing  the 
experiments . 
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CHAPTER  3 

EXPERIMENTAL  DETERMINATION  OF  RADIANT  ENERGY  TRANSFER 
FROM  A  LABORATORY  PLASMA 

\7hen  measuring  the  radiant  energy  transfer  from  a  plasma  gen¬ 
erated  in  a  plasma  torch,  that  portion  of  the  plasma  emitting  the 
radiation  being  measured  must  be  defined.  The  luminous  tongue  of 
gas  which  protrudes  from  the  exit  nozzle  of  the  plasma  torch  is  re¬ 
ferred  to  as  the  plasma  jet.  The  model  used  to  determine  the  total 
radiation  from  the  plasma  jet  is  described  in  Section  3»1*  The  me¬ 
thod  of  making  the  pointwise  radiation  measurements  along  the  axis 
of  the  plasma  jet  is  discussed  in  Section  3 • 5 • 

In  addition  to  the  radiation  emitted  from  the  plasma  jet,  ra¬ 
diation  is  emitted  from  the  plasma  contained  within  the  exit  nozzle 
of  the  plasma  torch,  which  is  referred  to  as  the  hot  core  of  plasma 
within  the  exit  nozzle.  The  measurements  used  to  determine  the  amount 
of  radiant  energy  which  is  actually  transferred  from  the  plasma  torch 
during  operation  are  described  in  Section  3»8»  This  radiant  energy 
transfer  is  referred  to  as  the  actiial  radiative  heat  transfer  from 
the  plasma  and  plasma  torch. 
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3.1  MODEL  FOR  TOTAL  RADIATION  FROM  A  PLASMA  JET 

The  luminous  portion  of  the  plasma  which  protrudes  from  the 
plasma  torch  is  referred  to  as  the  plasma  jet.  The  plasma  jets  studied 
in  these  experiments  were  roughly  conical  in  shape  resembling  flames. 
For  a  transparent  plasma,  the  intensity  of  the  radiation  emitted  from 
the  plasma  depends  on  the  volume  of  plasma  observed.  If  the  plasma 
jet  could  be  detached  from  the  plasma  torch  and  placed  at  the  center 
of  a  sphere,  the  radiation  emitted  from  the  plasma  jet  would  be 
distributed  uniformly  over  the  sphere.  Measurements  on  combustion 
flames  demonstrate  this  fact  (Eef.  3). 

In  order  to  represent  the  radiation  emitted  from  the  plasma 
jet,  the  plasma  jet  is  imagined  to  lie  at  the  center  of  a  sphere 
as  shown  in  Figure  3-1.  The  radius  of  this  sphere  is  taken  as  the 
distance  from  the  axis  of  the  plasma  jet  to  the  bolometer  sensing 
element.  The  total  radiation  emitted  from  the  plasma  jet  is  defined 
as  the  radiation  crossing  the  sphere,  assuming  that  this  radiation 
is  distributed  uniformly  over  the  sphere.  To  calculate  the  total 
radiation  from  the  plasma  jet,  the  radiative  flux  at  the  bolometer 
element  is  multiplied  by  the  area  of  the  sphere  surrounding  the 
plasma  jet. 
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Fig.  3-1  Model  Used  to  Calculate  Total  Radiation  from  Plasma 
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3.2  TOTAL  RADIATION  FROM  AN  ARGON  PLASMA  JET 

All  of  the  total  radiation  measurements  were  made  with  the  bolo¬ 
meter  and  plasma  torch  positioned  as  shown  in  Figure  3-2.  The  bolometer 
was  positioned  so  that  a  normal  to  the  surface  of  the  bolometer  sensing 
element  would  intersect  the  axis  of  the  plasma  jet  at  a  right  angle  at 
the  base  of  the  plasma  jet.  The  distance  from  the  bolometer  sensing 
element  to  the  axis  of  the  plasma  jet  was  ItO  3/8  in.  The  bolometer  was 
calibrated  in  this  position  as  discussed  in  Appendix  B.  In  all  of  the 
experiments  the  plasma  jet  was  exhausted  to  the  atmosphere. 

Radiation  readings  were  recorded  over  the  maximum  operating  range 
of  the  plasma  torch  by  varying  the  current  to  the  plasma  torch  incre¬ 
mentally.  The  current  was  usually  changed  in  50  amp  increments.  A 
minimum  of  one  minute  was  allowed  after  each  change  in  the  current  to 
permit  the  plasma  torch  and  all  instruments  to  reach  a  state  of  steady 
operation.  The  net  power  input  to  the  plasma  was  calculated  for  each 
radiation  reading.  The  total  radiation  from  the  plasma  jet  was  calcu¬ 
lated  as  discussed  in  Section  3.1.  Sample  calculations  are  given  in 
Appendix  C. 

The  results  of  the  measurements  of  the  total  radiation  from  an 
argon  plasma  jet  are  indicated  in  Figure  3-3  by  the  open  circles.  The 
total  radiation  from  the  plasma  jet  calculated  for  each  radiation  read¬ 
ing  is  expressed  as  a  percentage  of  the  net  power  input  to  the  plasma 
at  that  radiation  reading.  The  experimental  curve  has  been  drawn  by 
hand  as  the  best  apparent  representation  of  the  data. 

The  effect  of  the  flow  rate  of  argon  on  the  total  radiation  from 


16 


AEDC-TDR-64-233 


BOLOMETER  PLASMA  TORCH 

UNIT 


Fig.  3-2  Arrangement  for  Total  Radiation  Measurements 
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an  argon  plasma  jet  is  shown  in  Figure  3-4o  The  curves  are  similar 
in  shape  and  the  data  is  reproducible.  It  was  observed  that  for  an 
argon  flow  rate  of  2.44  Ibs/hr  the  plasma  jet  was  quiet  at  all  power 
levels  and  this  was  interpreted  to  be  laminar  flow.  Also  it  was  noted 
that  these  quiet,  laminar  jets  were  much  longer  than  the  turbulent 
jets  (higher  argon  flow  rates);  indicating  that  turbulent  mixing  between 
the  plasma  and  the  surrounding  air  has  a  quenching  effect.  At  inter¬ 
mediate  and  low  power  inputs  for  the  flow  rate  3.97  lbs/hr  the  jet  was 
laminar;  for  high  power  inputs  the  jet  was  turbulent.  Only  for  the 
lowest  power  inputs  was  the  jet  laminar  for  an  argon  flow  rate  of 
5.51  lbs/hr„  For  all  other  power  inputs  and  flow  rates  shown  in  Figure 
3-4  the  jet  was  turbulent. 

A  parameter  used  by  experimentalists  to  describe  the  state  of  the 
plasma  jet  is  the  average  enthalpy  input  per  pound  of  gas.  This 
quantity  is  defined  as  the  net  power  input  divided  by  the  gas  flow  rate, 
and  is  represented  by  the  45°  line  shown  in  Figure  3-4a.  This  line  is 
labeled  "radiation  loss  neglected".  To  take  into  account  the  radiation 
loss,  the  following  expression  was  used  to  calculate  the  enthalpy  rise? 

Net  Energy  Input  -  Radiant  Energy  Loss  =  Enthalpy  Rise.  All  the 
data  points  shown  in  Figure  3-4  were  used  to  calculate  the  enthalpy  rise 
and  the  results  are  seen  in  Figure  3-4a.  There  appear  to  be  two  distinct 
sets  of  data  points  -  one  for  the  laminar  plasma  jet  and  the  other  for 
the  turbulent  plasma  jet. 

In  defining  the  enthalpy  rise  many  quantities  have  been  neglected; 
such  as,  kinetic  energy  of  the  gases,  viscous  dissipation,  heat  conduction 
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etCo  It  is  believed  that  the  neglect  of  viscous  dissipation  may  account 
for  laminar  curve  being  above  the  turbulent  curve,,  Viscous  dissipation  is 
higher  for  turbulent  flows  than  for  corresponding  laminar  flows.  However 
it  is  not  claimed,  but  only  weakly  suggested  that  viscous  dissipation 
may  account  for  most  of  the  spread  between  the  laminar  and  turbulent 
curves  in  Figure  3-4a, 

3.3  TOTAL  RADIATION  FROM  A  NITROGEN  PLASMA  JET 

Measurements  of  the  total  radiation  from  a  nitrogen  plasma  jet  were 
made  in  the  same  manner  as  the  measurements  of  the  total  radiation  from 
an  argon  plasma  jet. 

The  operation  of  the  plasma  torch  with  nitrogen  was  much  more 
difficult  than  with  argon.  While  generating  a  nitrogen  plasma  the 
current  recorder  trace  was  unsteady,  making  it  possible  to  read  the 
value  of  the  current  only  to  within  +  30  amps.  Due  to  thh^  high  power 
level  necessary  to  generate  a  nitrogen  plasma,  the  nozzle  inserts  and 
"0"-ring  water  seals  were  short-lived.  The  nitrogen  plasma  jet  was 
much  different  in  appearance  than  the  argon  plasma  jet,  being  two  or 
three  times  longer  and  ragged  in  appearance.  The  nitrogen  plasma  jet 
was  slightly  yellow  in  color  as  opposed  to  the  bluish-shite  argon  plasma 
jet  and  was  not  nearly  as  intense. 

The  results  of  the  total  radiation  measurements  from  a  nitrogen 
plasma  jet  are  shown  in  Figure  3-5,  A  straight  line  has  been  used  to 


20 


2 


4 


10 


6  8 
NET  POWER  INPUT  (BTU/ 

Fig.  3-4  Effect  of  Flow  Rate  on  Toto 


Cooling  water  flew  rate  810  Ib/hr. 
Operating  voltage  range  19.3  -25. 9  volts 
Operating  current  ronge  ISO**  397  amps. 


□  Ascending  power  levels 
■  Descending  power  levels 

D.  Evans  —I 

Gas  Dynamics  Laboratory  q 

^  ^Northw®st®|n  University^  ^  ^  |  ^3 

12  14  16  18 

r.xlO"3) 


Radiation 


54-233 


ENTHALPY  RISE  CBTU/lbxlO3) 


AEDC-TDR-64-233 


.6  .8  1,0  1,2  1.4  1.6  IB  20  22  2.4  26  2B  30 

NET  ENERGY  INPUT  (BTU/lb  x  ICT3) 


Fig.  3-4a  Effect  of  Radiation  on  Enthalpy  Rise 


22 


PERCENT  NET  POWER  RADIATED 


AEDC-TDR-64-233 


NET  POWER  INPUT  TO  GAS  (Btu/hr  x  10"*) 

Fig.  3-5  Total  Radiation  from  a  Nitrogen  Plasma 
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represent  the  data  due  to  the  expanded  scale  along  the  ordinate  and 
the  increased  error  in  the  power  measurements. 

3.4  TOTAL  RADIATION  FROM  AN  ARGON  PLASMA  JET  SEEDED  WITH  CARBON 
PARTICLES 

An  argon  plasma  jet  was  seeded  with  finely  divided  carbon  parti¬ 
cles  in  an  attempt  to  increase  the  optical  thickness  of  the  plasma  jet. 
Carbon  black  manufactured  by  Cabot  Corporation  under  the  brand  name 

Supercarbovar  was  used  as  the  seeding  material.  The  carbon  particles 

8 

used  have  a  mean  diameter  of  I.I4  x  10“  m.  (Ref.  U).  A  Sylco  CCC  Mark 
VU  fine  powder  feeder  manufactured  by  Sylvester  and  Company  was  used 
to  feed  the  carbon  to  the  plasma  jet. 

The  total  radiation  from  the  seeded  argon  plasma  jet  was  measured 
with  carbon  being  fed  to  the  plasma  jet  at  four  different  mass  flow 
rates.  Total  radiation  measurements  with  no  carbon  being  fed  to  the 
plasma  jet  were  taken  preceding  each  reading  with  carbon  particles  pre¬ 
sent  in  the  jet.  In  this  manner  eight  readings  were  taken  at  each  of 
six  different  power  levels.  The  measurements  ware  made  in  the  same 
manner  as  the  total  radiation  measurements  from  the  argon  and  nitrogen 
plasma  jets. 

The  results  of  these  measurements  are  shown  in  Figure  3-6.  Fluctua¬ 
tions  in  the  voltage  across  the  plasma  torch  and  a  slight  amount  of 
electrical  noise  in  the  current  recorder  trace  due  to  the  seeder  account 
for  the  fact  that  each  group  of  data  points  is  not  at  precisely  the  same 
net  power  input  level.  The  plasma  jet  took  on  a  slight  yellow  coloration 
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Fig.  3-6  Total  Radiation  from  an  Argon  Plasma  Seeded  with  Carbon  Particles 
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with  carbon  particles  present,  an  effect  which  was  more  pronounced  at 
the  higher  carbon  flow  rates.  Carbon  was  observed  leaving  the  tip  of 
the  plasma  jet  as  a  fine  black  smoke,  giving  the  jet  an  appearance 
similar  to  a  sooty  flame. 

3.5  POINTWISE  RADIATION  MEASUREMENTS 

Measurements  of  the  radiative  flux  at  the  surface  of  the  plasma 
jet  were  made  by  projecting  an  image  of  the  plasma  jet  in  the  plane  of 
the  bolometer  sensing  element.  These  measurements  are  referred  to  as 
pointwise  radiation  measurements. 

The  optical  system  shown  in  Figure  3-7  was  used  to  project  an 
image  of  the  plasma  jet  in  the  plane  of  the  bolometer  element.  Two¬ 
fold  magnification  was  used  for  convenience  to  provide  a  larger  image 
with  which  to  work  and  to  increase  the  resolution  of  the  bolometer. 

The  mirror  used  is  a  10  in.  diameter  front-surfaced  paraboloid  having  a 
in.  focal  length.  The  optical  system  was  aligned  using  a  transparent 
focusing  grid  and  a  light  source  as  an  object.  The  calibration  of  the 
bolometer  for  use  with  this  optical  system  is  discussed  in  Appendix  B. 

3.6  DISTRIBUTION  OF  RADIATIVE  FLUX  ALONG  THE  AXIS  OF  AN  ARGON  PLASMA  JET 

With  the  arrangement  shown  in  Figure  3-7,  the  bolometer  was  moved 
horizontally  across  the  image  of  the  plasma  jet  until  a  maximum  output 
signal  was  obtained.  This  position  corresponds  to  the  location  of  the 
axis  of  the  plasma  jet  in  its  image.  Maintaining  the  current  to  the 
plasma  torch  at  a  constant  value,  the  image  was  traversed  vertically 

/ 
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along  the  axis*  The  bolometer  output  signal  and  the  voltage  drop 
across  the  plasma  torch  were  recorded  simultaneously  at  l/O  in.  incre¬ 
ments.  Ascending  and  descending  traverses  were  taken  at  each  of  four 
power  levels.  The  data  from  four  of  these  eight  traverses  are  shown 
in  Figure  3-8.  The  net  power  input  to  the  plasma  was  calculated  for 
each  radiation  reading.  These  values  were  then  averaged  to  compute  the 
average  net  power  for  each  traverse®  The  standard  deviation  from  this 
average  was  also  calculated,  and  was  less  than  one  per  cent  in  all  but 
one  case.  In  Figure  3-8,  the  argon  flow  rate  at  the  three  higher  power 
levels  was  8.5  Ibs/hr;  the  flow  rate  at  the  lowest  power  level  was  9.5 
lbs/hr.  Figure  3-8a  shows  this  data  plotted  logarithmically . 

If  the  maximum  output  signal  along  the  axis  of  the  image  of  the 
plasma  jet  is  taken  as  the  base  of  the  plasma  jet,  the  distribution  of 
the  radiative  flux  along  the  axis  of  the  plasma  jet  is  represented 
closely  by  the  expression 


F  =  F 


max 


where  F  is  the  radiative  flux  at  any  point  x  along  the  axis  of  the 

2 

plasma  jet  (Btu/hr  ft  ) 

F  is  the  maximum  radiative  flux,  i,e.«  the  radiative  flux  at  the 
max  ’  s 

2 

base  of  the  plasma  jet  (Btu/hr  ft'  ) 
x  is  the  distance  from  the  tip  of  the  plasma  jet  to  the  observa¬ 
tion  point  (in) 

I.  is  the  length  of  the  plasma  jet  from  base  to  tip  (in) 
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Fig.  3-8  Pointwise  Radiation  Measurements  of  Argon  Plasma 
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Both  the  length  of  the  plasma  jet  and  the  maximum  radiative  flux  at  the 
base  of  the  plasma  jet  depend  on  the  net  power  input  to  the  plasma  and 
the  flow  rate  of  the  gas.  To  convert  the  bolometer  output  signals  given 
in  Figure  3-8  to  values  of  radiative  flux  at  the  surface  of  the  plasma 
jet,  the  calibration  value  of  l.S?l*  x  10^  Btu/hr  ft2  volt  can  be  used 
(see  Appendix  B). 

3.7  CALCULATION  OF  TOTAL  RADIATION  USING  POINTWISE  RADIATION  DATA 

Tine  total  radiation  from  an  argon  plasma  jet  was  calculated  using 
the  pointwise  radiation  data  shown  in  Figure  3-8.  The  calculations  are 
based  on  the  model  shown  in  Figure  3-9.  The  plasma  jet  is  assumed  to  be 
conical  in  shape  having  a  base  with  diameter  equal  to  the  diameter  of 
the  exit  nozzle  opening,  9/l6  in.,  and  having  a  height  as  determined 
from  the  pointwise  radiation  measurements.  The  cone  is  imagined  to  be 
sectioned  into  frustums  of  equal  altitude,  such  that  each  pointwise 
radiation  measurement  lies  at  the  vertical  midpoint  of  its  corresponding 
frustum.  The  total  radiation  from  the  argon  plasma  jet  was  calculated 
by  multiplying  the  radiative  flux  at  the  vertical  midpoint  of  each  frus¬ 
tum  by  the  area  of  that  frustum,  and  summing  the  contributions  from  all 
frustums  of  the  cone.  Examples  of  these  calculations  are  given  in 
Appendix  C.  The  results  of  these  calculations  are  plotted  on  Figure 
3-3  as  the  closed  circles. 

3.8  ANGULAR  DISTRIBUTION  OF  RADIATION  FROM  PLASMA  TORCH 

The  angular  distribution  of  the  radiation  from  the  plasma  torch 
was  determined  by  mounting  the  plasma  torch  on  a  small  rotary  machinist's 
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Fig.  3-9  Model  used  to  Calculate  Total  Radiation  from  Plasma  using  Pointwise  Data 
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table,  as  shown  in  Figure  3-10,  The  plasma  torch  was  positioned  on  the 
table  so  that  the  center  of  rotation  was  on  the  axis  of  the  plasma  jet 
at  its  base.  The  bolometer  sensing  element  was  located  U7  5/8  in.  from 
the  center  of  rotation  in  the  plane  of  rotation. 

The  experiment  was  performed  with  the  plasma  torch  generating  an 
argon  plasma.  The  angular  distribution  of  the  radiation  from  the  plasma 
torch  was  determined  by  rotating  the  plasma  torch  in  angular  increments. 
Radiation  readings  were  recorded  at  5°  increments  at  the  smaller  viewing 
angles  and  at  10°  increments  at  the  larger  viewing  angles.  The  voltage 
drop  across  the  plasma  torch  was  recorded  simultaneously  with  each 
radiation  reading.  The  readings  were  begun  at  a  viewing  angle  of  0°, 
with  the  bolometer  sensing  element  lying  on  the  extended  axis  of  the 
plasma  jet..  The  plasma  torch  was  rotated  until  the  bolometer  output 
signal  fell  to  the  noise  level  in  the  shadow  of  the  plasma  torch.  The 
air  current  from  a  household  fan  was  directed  transverse  to  the  line 
from  the  bolometer  sensing  element  to  the  center  of  rotation  of  the 
plasma  torch  to  prevent  convective  heating  at  the  bolometer. 

The  angular  distribution  of  the  radiation  was  determined  at  four 
power  levels.  The  results  of  these  measurements  are  shown  in  Figure 
3-11.  The  net  power  input  to  the  plasma  was  calculated  at  each  radia¬ 
tion  reading  and  these  values  were  averaged  over  each  run  to  compute  the 
average  net  power. 

3.9  ACTUAL  RADIATIVE  HEAT  TRANSFER  FROM  AN  ARGON  PLASMA  AND  THE  PLASMA 
TORCH 

With  the  plasma  torch  operating  at  a  given  set  of  conditions,  there 
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Arrangement  for  Measurement  of  Angular  Distribution  of  Radiation 
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Fig.  3-11  Angular  Distribution  of  Radiation  from  Plasma  Torch 


A  EDC-T  DR-64-233 


is  a  definite  amount  of  radiant  energy  transferred  to  the  surroundings. 
This  is  referred  to  as  the  actual  radiative  heat  transfer  from  the 
plasma  and  plasma  torch. 

The  data  of  the  angular  distribution  of  the  radiation  from  the 
plasma  torch  shown  in  Figure  3-11  were  used  to  calculate  the  actual 
radiative  heat  transfer  from  an  argon  plasma  and  the  plasma  torch. 

These  calculations  are  based  on  the  model  shown  in  Figure  3-12.  The 
bolometer  is  imagined  to  lie  on  the  surface  of  a  sphere  whose  center  is 
the  center  of  rotation  of  the  plasma  torch.  At  a  viewing  angle  of  0° 
the  bolometer  sensing  element  lies  on  the  south  pole  of  this  sphere. 

The  surface  is  imagined  to  be  divided  into  ring-like  segments,  such 
that  each  radiation  reading  lies  in  the  latitudinal  center  of  its 
corresponding  segment.  The  radiation  crossing  each  segment  of  the 
surface  was  calculated  by  multiplying  the  flux  at  that  segment  by  the 
area  of  the  segment.  The  actual  radiative  heat  transfer  was  taken  as 
the  sum  of  the  radiation  crossing  all  segments  of  the  spherical  surface. 
Examples  of  these  calculations  are  shown  in  Appendix  C. 

The  results  of  these  calculations  are  shown  in  Figure  3-13  •  The 
actual  radiative  heat  transfer  from  the  argon  plasma  and  the  plasma 
torch  is  expressed  as  a  percentage  of  the  average  net  power  input  to 
the  plasma  for  each  determination. 

3.10  DISCUSSION  OF  RESULTS 

The  efficiency  of  the  plasma  torch  is  defined  as  the  ratio  of  the 
net  power  input  to  the  plasma  to  the  gross  electrical  power  input  to 
the  plasma  torch.  While  generating  an  argon  plasma,  the  efficiency  of 
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Fig.  3-12  Model  for  Calculation  of  Actual  Radiative  Heat  Transfer 
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Fig.  3-13  Actual  Radiative  Heat  Transfer  from  Argon  Plasma  and  Plasma  Torch 
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the  plasma  torch  remained  nearly  constant  at  0.5.  For  the  forty  data 
points  recorded  in  the  determination  of  the  total  radiation  from  an 
argon  plasma  (Figure  3-3)  s  the  efficiency  varied  from  0.1j8  to  0o$k  with 
an  average  of  0.5>1®  The  efficiency  of  the  plasma  torch  while  generating 
an  argon  plasma  is  taken  as  constant  at  0 .5  over  the  range  of  power 
levels  investigated. 

The  primary  source  of  errors  in  the  experiments  described  above 
is  the  power  measurements.  The  cooling  water  flow  rate  was  observed  to 
fluctuate  in  a  random  manner  approximately  ±  5  per  cent.  The  dead  band 
of  the  current  recorder  is  1  per  cent.  The  accuracy  of  the  power  volt¬ 
meter  is  2  per  cent.  This  gives  an  overall  possible  error  of  ±  11  per 
cent  in  the  determination  of  the'  net  power  input. 

Referring  to  Figure  3-3 »  most  of  the  scatter  of  the  data  can  be 
attributed  to  the  errors  incurred  in  the  power  measurements.  There  is, 
however,  an  additional  source  of  error  worth  considering.  As  mentioned 
in  Section  1.3,  the  use  of  the  net  power  input  to  the  plasma  as  the 
independent  variable  requires  the  implicit  assumption  that  the  net  power 
input  to  the  plasma  uniquely  determines  the  temperature  distribution  in 
the  plasma  jet.  The  fact  that  the  radiation  emitted  from  the  plasma  jet 
shows  a  strong  correlation  with  the  net  power  input  to  the  plasma  in¬ 
directly  verifies  this  assumption. 

The  total  radiation  from  an  atmospheric  argon  plasma  jet  is  signi¬ 
ficant,  vaiying  between  1  per  cent  and  18  per  cent  of  the  net  power  input 
to  the  plasma  over  the  range  of  power  levels  Investigated®  At  a  constant 
net  power  input,  the  total  radiation  from  an  argon  plasma  jet  increases 
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with  decreasing  flow  rate,  due  to  the  higher  temperatures  generated  in 
the  plasma  jet  at  the  lower  flow  rates. 

The  total  radiation  from  an  atmospheric  nitrogen  plasma  jet  is 
less  than  1  per  cent  of  the  net  power  input  to  the  plasma  over  the  range 
of  power  levels  investigated.  Since  nitrogen  is  a  diatomic  gas,  most 
of  the  energy  input  to  the  gas  is  expended  in  dissociating  the  molecules. 
As  a  result,  the  population  of  the  free  charged  particles,  which  is  the 
source  of  the  radiation  emitted  from  a  plasma,  is  relatively  small.  For 
this  reason,  the  total  radiation  from  a  nitrogen  plasma  jet  is  at  least 
an  order  of  magnitude  smaller  than  the  total  radiation  from  an  argon 
plasma  jet  at  the  same  net  power  input  to  the  plasma.  The  total  radia¬ 
tion  from  an  atmospheric  nitrogen  plasma  jet  may  be  considered  insigni¬ 
ficant  for  most  purposes  over  the  range  of  power  levels  investigated. 

The  calculation  of  the  total  radiation  from  the  argon  plasma  jet 
using  the  pointwise  radiation  data  agrees  reasonably  well  with  the 
direct  measurements.  The  fact  that  this  independent  method  of  determin¬ 
ing  the  total  radiation  from  an  argon  plasma  jet  agrees  with  the  results 
of  the  direct  measurements  tends  to  confirm  the  results  of  the  direct 
measurements  of  the  total  radiation  from  an  argon  plasma  jet. 

The  question  of  the  effect  of  the  presence  of  finely  divided  car¬ 
bon  particles  on  the  total  radiation  from  an  argon  plasma  jet  must  be 
left  unanswered.  The  results  of  the  experiments  reported  here  reveal 
no . measurable  effect  due  to  the  carbon.  Further  investigations  in  this 
area  were  not  attempted. 

The  data  on  the  angular  distribution  of  the  radiation  from  the 
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plasraa  and  plasma  torch  reveal  the  presence  of  a  hot  core  of  plasma 
within  the  exit  nozzle  which  radiates  intensely.  As  the  viewing  angle 
decreases,  more  of  the  radiation  from  this  hot  core  of  plasma  is  detected 
by  the  bolometer.  Thus,  the  intensity  of  radiation  increases  with 
decreasing  viewing  angle.  The  attenuation  of  the  radiation  as  the  view¬ 
ing  angle  decreases  from  10°  to  0°  is  evidence  of  absorption  of  the  radia¬ 
tion  from  the  hot  core  of  plasma  by  the  plasma  jet.  At  a  viewing  angle 
of  0°,  the  radiation  emanating  from  within  the  exit  nozzle  must  traverse 
the  entire  length  of  the  plasma  jet  before  being  detected  by  the  bolo¬ 
meter.  Thus,  the  argon  plasma  generated  in  these  experiments  begins  to 
exhibit  a  measureable  absorptivity  at  a  depth  of  about  one  inch. 

The  determination  of  the  actual  radiative  heat  transfer  from  the 
argon  plasma  and  the  plasma  torch  indicates  that  radiation  is  a  signi¬ 
ficant  mechanism  of  heat  transfer  from  an  atmospheric  argon  plasma 
generated  in  a  plasma  torch.  Assuming  a  constant  efficiency  of  0.5  for 
the  plasma  torch,  these  data  can  be  used  to  calculate  an  overall  energy 
balance  about  the  plasma  torch  including  radiative  losses  from  the  argon 
plasma. 
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CHAPTER  4 

THEORETICAL  PREDICTION  OF  RADIATION  EMITTED  FROM  A 
PARTIALLY  IONIZED  GAS 

4.1  SOURCE  OF  CONTINUUM  RADIATION 

The  major  portion  (energetically)  of  the  radiation  emitted  "by  a 
hot,  partially  ionized  gas  (plasma)  at  atmospheric  pressure  is  a  con¬ 
tinuous  function  of  frequency.  The  source  of  this  continuum  radiation 
is  the  motion  of  free  charges  in  the  plasma.  It  is  known  from  electro¬ 
dynamics  that  a  free  charge,  when  accelerated,  emits  electromagnetic 
radiation.  In  a  plasma  the  free  charges  (electrons)  are  accelerated 
as  they  collide  with  other  particles  (atoms,  ions,  and  electrons), 
thus  emitting  electromagnetic  radiation. 

According  to  quantum  mechanics,  the  energy  of  electromagnetic 
radiation  is  proportional  to  frequency  through  Planck's  constant 


E  *  h-v 


(4-1) 


where  h  is  Planck's  constant,  6.63  x  10~^  erg  sec 
■V  is  the  frequency  (sec“l) 

E  is  the  energy  (ergs) 

Since  energy  is  conserved,  an  electron  undergoing  a  collision  has 
some  of  its  kinetic  energy  converted  to  radiation.  The  frequency  of 
the  radiation  emitted  by  the  electron  is  proportional  to  the  energy  loss 
of  the  electron.  A  free  electron,  upon  colliding  with  another  particle, 
can  undergo  a  transition  to  any  lower  energy  state;  thus  the  radiation 
emitted  from  such  collisions  is  continuous  in  frequency.  This  continuous 
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radiation  is  sometimes  referred  to  as  bremsstrahlung. 

4.2  FREE-FREE  TRANSITIONS 

An  electron,  upon  colliding  with  an  ion,  lias  some  of  its  kinetic 
energy  converted  to  radiation,,  If  the  energy  loss  of  the  electron  is 
not  too  great,  it  leaves  the  influence  of  the  ion's  Coulombic  field  and 
maintains  its  status  as  a  free  electron*  This  is  an  example  of  a  colli¬ 
sion  classified  as  a  free-free  transition. 

The  radiation  emitted  by  an  electron  undergoing  a  free-free  transi¬ 
tion  in  a  collision  with  an  ion  has  been  analyzed  by  Oster  (Ref.  5). 

His  treatment  is  a  classical  one,  in  which  the  radiant  energy  loss  of 
the  electron  is  ignored  in  solving  the  mechanics  of  the  collision*  As 
a  result  of  this  assumption,  the  electron  describes  a  hyperbolic  orbit 
in  a  reference  frame  attached  to  the  ion*  A  Fourier  analysis  of  the 
components  of  the  electronic  acceleration  yields  the  classical  emission 
coefficient  of  a  free -free  electron-ion  collision* 

Because  of  the  mathematical  complexity  of  the  expression  thus  ob¬ 
tained,  one  of  two  approximations  to  the  hyperbolic  orbit  of  the  elec¬ 
tron  is  commonly  made.  In  one  case,  the  hyperbolic  orbit  is  approximated 
by  a  straight  line,  in  the  other  by  a  parabola.  Each  of  these  approxi¬ 
mations  yields  results  which  are  valid  over  certain  frequency  ranges. 
Using  Kramers'  parabolic  approximation  to  the  electronic  orbit  (Ref.  6), 
a  simple  expression  for  the  emission  coefficient  is  obtained  from  a 
Fourier  analysis  of  the  electron's  acceleration  (Refs.  7,  8) 


43 


AEDC-TDR-64-233 
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a-2) 


where  e  is  the  electronic  charge,  U.8l  x  lO”'^  statcoul 

*»28 

m  is  the  electronic  mass,  9*12  x  10“  gms 
c  is  the  velocity  of  light,  3  x  10^  cm/sec 
Z  is  the  integral  ionic  charge 

nQ,  n^  are  the  electron,  ion  number  densities  respectively  (cm  ) 
T  is  the  temperature  (°K) 

—27 

h  is  Planck's  constant,  6.63  x  10  erg  sec 
k  is  Boltzmann's  constant,  1*38  x  10 erg/°IC 
v  is  the  frequency  of  radiation  (sec-’1') 

cy  is  the  radiant  energy  of  frequency  v  emitted  in  all  directions 
per  unit  volume  of  gas  (ergs/cm  J) 

Evaluating  the  constant  term,  this  equation  reduces  to 


=  6.8  X  10-39  Z2 


04-3) 


Integrating  this  expression  over  all  frequencies  gives  the  power  radiated 
per  unit  volume  of  gas  at  temperature  T 


P(T) 


& 

/«. 

0 


dv 


(Wi) 


where  P(T)  is  the  power  radiated  per  unit  volume  (ergs/ cm'1  sec).  Substi¬ 
tuting  Eq.  L|  —  3  into  Eq.  I1-I4  yields 


44 


AEDC-TDR- 64-233 


P(T)  a  J 5.8  x  |0“88  Z2 

5  T,,z 

Since  the  ionic  charge,  number  densities, 
of  frequency,  this  gives 


e-^/kT  dv  (1,-5) 

and  temperature  are  independent 


run.-  /.oo  ,  „  _ 

P(T)  =  6.8  x  IO~ 38  Z2  w  /  firhi//kT  ^ 

o 

=  6.8  X  I0“38  Z2  (~  )  ff-hv/kTl® 

1  n  *0 

=  6.8  x  10“ 38  £  Z2  nenj  T,/2 

o 

Evaluating  the  constant  term  and  expressing  the  result  in  watts/cnr  gives 
(Ref.  9) 

P(T)  =  1.4  XIO"34  Z2  nenj  Tl/2  (1,-6) 

In  applying  the  Kramers  approximation  to  the  work  discussed  in 
Chapter  5,  the  limits  of  integration  used  were  the  upper  and  lower  cut¬ 
off  frequencies  of  the  sodium  chloride  window.  Substituting  these  limits 
into  Eq*  1,-5  gives 


P (T )  =  1.4  xlO-34  Z2nen,-  T,/2  [ff'84o/T  -  ^-sa.eoo/T J  (i,.7) 

This  result  is  plotted  as  a  function  of  temperature  in  Figure  5-1 9. 

The  values  of  electron  and  ion  number  densities  were  taken  from  the 
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tables  computed  by  Drellishak,  et  al.  (Ref.  10)  for  the  equilibrium 
composition  of  atmospheric  argon  plasma.  Calculations  were  performed 
at  thousand  degree  intervals  from  10,000°K  to  23jOOO°K.  Two  calcula¬ 
tions  were  perfomed  at  each  temperature  ,  one  for  the  radiation  from 
collisions  with  the  first  ions  (Z  =  l)  and  one  for  the  radiation  from 
collisions  with  the  second  ions  (Z  =  2).  These  two  contributions  were 
added.  The  effect  of  radiation  from  electrons  colliding  with  second 
ions  calculated  from  Eq.  Jh—7  is  negligible  below  18,000°K*  Sample 
calculations  are  shown  in  Appendix  D. 

4.3  QUANTUM  CORRECTIONS  TO  CLASSICAL  RESULTS 

The  basic  assumption  which  allows  a  classical  treatment  of  the 
radiation  emitted  by  an  electron  upon  colliding  with  an  ion  is  the  dis¬ 
regard  for  the  energy  lost  by  the  electron  to  radiation  during  the  en¬ 
counter.  A  comprehensive  quantum  mechanical  analysis  of  the  electron- 
ion  encounters  which  takes  account  of  the  energy  loss  of  the  electron 
was  achieved  by  Sommerfeld  in  his  book  (1939)  on  atomic  spectra.  Greene 
(Ref.  11)  gives  a  correction  to  Kramers*  classical  equation  in  the  form 
of  a  multiplicative  factor  which  adjusts  the  classical  Kramers  equation 
to  agree  with  Sommerfeld  *  s  analysis 

P(T)  =  1.4  XIO"34  Z2  ne rij  Tl/2  g(T)  Oi-Q) 

The  multiplicative  factor,  g(T),  is  a  frequency  averaged  Gaunt  factor 
and  is  a  slowly  varying  function  of  temperature.  A  plot  of  Eq.  l|-8, 
calculated  using  values  of  g(T)  given  by  Greene  is  shown  on  Figure  5-19 
as  the  Sommerfeld  solution. 
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4.4  FREE-BOUND  TRANSITIONS 

Consider  now  the  radiation  emitted  by  a  free  electron  as  it 
collides  with  and  is  captured  by  an  ion.  Such  encounters  are  called 
free-bound  transitions.  A  free -bound  transition  occurs  when  a  free 
electron  loses  a  sufficient  amount  of  its  energy  and  is  captured  by  an 
ion.  The  radiation  emitted  by  electrons  undergoing  free-bound  transi¬ 
tions  is  often  referred  to  as  recombination  radiation. 

One  of  the  basic  tenets  of  quantum  theory  is  that  electrons  are 
bound  to  their  nuclei  in  discreet  energy  levels  or  orbits.  Elementary 
quantum  theory  predicts  that  the  energy  difference  between  successively 
more  energetic  orbits  decreases,  approaching  zero  as  the  orbital  energy 
approaches  the  ionization  energy. 

A  free  electron  which  undergoes  a  free -bound  transition  becomes 
bound  to  the  ion  in  one  of  these  discreet  energy  levels  or  orbits.  A 
calculation  of  the  radiation  emitted  by  electrons  experiencing  free- 
bound  transitions  would  involve  a  sum  over  these  discreet  levels  or 
states. 

In  a  hot  plasm,  however,  the  ions  are  excited  so  that  the  proba¬ 
bility  is  greatest  that  an  electron  will  become  bound  to  an  ion  in  one 
of  the  highest  energy  levels.  Olsen  (Ref.  12)  discusses  an  approxima¬ 
tion  which  Unsold  has  made  to  calculate  the  radiation  emitted  during 
free -bound  transitions. 

UnsOld  replaces  the  sum  over  discreet  states  by  an  integral. 

This  amounts  to  assuming  that  the  upper  orbits  of  an  ion  which  are 
available  to  accept  an  incoming  electron  are  so  energetically  dense 
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that  the  substitution  of  an  integral  is  valid,  Unsttld  adds  the  radia- 
tion  from  free-bound  transitions  to  that  from  free-free  transitions  as 
calculated  from  Kramers®  approximation  and  arrives  at  a  spectrum  for 
the  continuum  which  is  independent  of  frequency  (Ref,  12) 

-  6.8  x  10“ 45  Z2  (h-9) 

where  ev  is  the  radiant  energy  of  frequency  v  emitted  in  all  direc- 
tions  per  unit  volume  of  gas  (watt  sec/cnr) 

Z  is  the  integral  ionic  charge 

n  ,  n.  are  the  number  densities  of  the  electrons  and  ions, 

S  1 

respectively  (cm  ) 

T  is  the  temperature  (°K) 

Unsdld  Js  approximation ,  Eq„  h-9s  is  valid  up  to  a  cutoff  frequency, 
t/g  9  This  cutoff  frequency  is  defined  by  the  energy  difference  between 
the  effective  ionization  potential  and  the  energy  level  below  which  the 
assumption  of  dense  states  becomes  invalid.  The  theory  is  based  on 
recombination  into  hydrogen-like  states. 

The  total  power  radiated  from  the  gas  over  the  range  of  validity 
of  Unsdld  *s  approximation  is  easily  calculated  by  integrating  the 
emission  coefficient,  Eq,  li-9,  over  the  appropriate  frequency  interval 
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P(T)  =  f9  ev  dv 
© 


P(T)  =  6.8  x  I0"4§  Z2  neni 

T1'2  9  (1,-10) 

This  expression  gives  the  radiation  from  free-bound  and  free-free 
transitions  in  the  frequency  range  up  to  the  cutoff  frequency. 

Olsen  (Ref.  12)  has  attempted  to  measure  the  cutoff  frequency,, 
va  ,  above  which  the  UnsSld  approximation  is  invalid.  He  made  measure- 
ments  of  the  intensity  of  the  continuum  at  l6,000°K  in  an  atmospheric 
argon  plasma  as  a  function  of  frequency.  According  to  Unsold  *  s  approxi¬ 
mation,  a  plot  of  the  intensity  of  the  continuum  radiation  as  a  function 
of  frequency  should  be  flat  out  to  the  cutoff  frequency,  Vg  ,  and  then, 
according  to  Kramers’  approximation,  should  fall  off  exponentially  (Ref. 
8).  Olsen  apparently  never  reaches  the  cutoff  frequency  in  his  measure¬ 
ments,  although  it  appears  from  his  experimental  data  that  he  may  be 
only  slightly  above  it.  He  concludes  only  that  the  cutoff  frequency  in 
an  atmospheric  argon  plasma  at  16,000  K  is  less  than  3  x  10i4  sec”  . 

4.5  KRAMERS-UNSOLD  MODEL  FOR  CONTINUUM  RADIATION 

It  is  now  possible  to  formulate  a  simple  model  of  the  continuum 
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radiation  frora  a  plasraa  which  will  account  for  the  radiation  from  free- 
free  and  free -bound  transitions,,  This  is  done  by  matching  the  UnsBld 
approximation,  Eqc  li -9,  to  the  Kramers  approximation,  Eq.  h-3g  at  the 

cutoff  frequency,  and  then  integrating  both  approximations  over  the 

; 

appropriate  frequency  intervals. 

Details  of  the  calculations  are  shown  in  Appendix  D.  The  cutoff 

ll’  “=*1 

frequency  was  taken  as  constant  at  3  x  10  1  sec"  .  The  contribution  to 
the  continuum  radiation  below  the  cutoff  frequency  by  free-bound  and 
free-free  transitions  was  calculated  from  the  Unsold  approximation  as 
indicated  in  Eq.  Ij-10.  However,  in  order  to  match  the  experimental  work, 

1  j  -| 

the  lower  limit  of  integration  was  taken  as  0.175  x  10  see"  ,  the 
lower  limit  of  transmissibility  of  the  sodium  chloride  window.  The 
emission  coefficient  as  given  by  the  Kramers  approximation,  Eq.  5-3, 
was  forced  to  match  the  emission  coefficient  as  given  by  the  Unsdld 
approximation,  Eq.  5-9,  at  the  cutoff  frequency,  Vg  ,  at  each  tempera¬ 
ture.  This  was  done  by  adjusting  the  constant  term  in  the  Kramers 
approximation,  Eq.  5-3.  The  resulting  expression,  with  the  new  constant, 

was  then  integrated,  using  as  the  lower  limit  of  integration  the  cutoff 

lli 

frequency,  Vg  ,  and  as  the  upper  limit  11.2  x  10  sec  '  ,  which  is  the 
upper  limit  of  transmissibility  of  the  sodium  chloride  window.  These 
calculations  were  performed  at  thousand  degree  increments  from  10,000°K 
to  23,000°K,  to  cover  the  range  of  experimental  observation.  The  re¬ 
sulting  curve  is  shown  on  Figure  5-1 9  as  the  Kramers -Uns bid  approxima¬ 
tion. 


50 


CHAPTER  5 


AEDC-TDR-64-233 


RADIATION  FROM  AN  ATMOSPHERIC  ARGON  PLASMA  AS  A 
FUNCTION  OF  TEMPERATURE 


5.1  LOCAL  THERMODYNAMIC  EQUILIBRIUM 


The  radiation  emitted  from  a  plasma  will  he  a  function  of  the  ther 
modynamic  state  of  the  plasma  provided  the  plasma  is  in  local  thermody¬ 
namic  equilibrium.  Local  thermodynamic  equilibrium  implies  that  the 
various  components  of  the  plasma  are  at  the  same  temperature  at  each 
point  within  the  plasma. 

Knopp  (Ref.  1J>)  has  investigated  the  existence  of  local  thermo¬ 
dynamic  equilibrium  in  an  argon  plasma  jet  at  atmospheric  pressure. 

This  investigation  is  based  on  the  elastic  nature  of  collisions  between 
electrons  and  heavy  particles  in  the  plasma.  Compton  (Refs.  14,  15) 
has  shown  that  these  collisions  are  elastic  provided 


where  E  is  the  electric  field  intensity  (volts/cm) 
p  is  the  pressure  of  the  plasma  (mm  Hg) 

The  electrode-nozzle  spacing  used  in  these  experiments  was  about  0.25 
cm.  Using  a  representative  voltage  drop  of  25  volts,  the  criterion  of 
elastic  collisions  is  satisfied  for  an  atmospheric  argon  plasma  jet. 

The  argon  plasma  jet  may  be  considered  to  be  in  local  thermodynamic 
equilibrium  under  these  conditions  (Ref.  15). 

McGregor  (Ref.  15a)  states  that  this  criterion  may  be  inadequate 
for  flowing  systems  at  high  densities  such  as  are  encountered  in  a 
plasma  torch  exhausting  to  atmospheric  pressure.  He  mentions  that  a 
criterion  for  equilibrium  in  such  a  flow  system  should  at  least  consider 
the  dwell  time  of  the  particles  in  the  arc.  A  detailed  investigation 
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of  the  flow  between  the  electrodes  and  the  structure  of  the  arc  in  the 
plasma  torch  vised  in  this  report  was  not  undertaken.  The  plasmas 
studied  in  this  report  areassumed  to  he  in  local  thermodynamic  equili¬ 
brium. 

5.2  DETERMINATION  OF  TEMPERATURE  DISTRIBUTION  IN  AN  ARGON  PLASMA  JET 

The  temperature  distribution  in  the  argon  plasma  jet  was  deter¬ 
mined  by  a  method  previously  employed  by  several  investigators  (Refs. 

13,  16,  17,  18).  The  method  is  based  on  a  knowledge  of  the  intensity 
of  a  spectral  line  of  argon  as  a  function  of  temperature.  The  theoreti¬ 
cal  basis  for  this  method  of  temperature  determination  is  discussed  in 
Appendix  E. 

A  Perkin-Elmer  Model  12C  spectrometer  was  used  to  measure  the 
relative  intensity  of  the  radiation  from  the  plasma  jet  in  appropriate 
frequency  intervals.  The  arrangement  of  the  apparatus  used  ■'in  these 
measurements  is  shown  in  Figure  5-1.  An  f/4.5  optical  system  consist¬ 
ing  of  two  compound  lenses  mounted  on  an  optical  bench  was  used  to  pro¬ 
ject  an  image  of  the  plasma  jet  without  magnification  in  the  plane  of 
the  entrance  slit  of  the  spectrometer.  The  plasma  torch  was  mounted 
vertically  on  a  small  rotary  machinist's  table  with  one  lead  screw  of 
the  table  oriented  perpendicular  to  the  axis  of  the  optical  system. 

This  lead  screw  was  driven  by  a  reversible  1  r.p.m.  synchronous  motor. 
With  this  mechanism,  the  plasma  jet  was  traversed  past  the  spectroscope 
at  0.100  in/min.  Two  razor  blades  spaced  0,005  in.  apart  were  used 
to  provide  a  horizontal  entrance  slit  directly  in  front  of  the  vertical 
entrance  slit  of  the  spectrometer.  This  combination  of  entrance  slits 
allowed  the  spectroscope  to  view  only  a  small  area  of  the  image  of  the 
plasma  jet. 

The  traverses  were  taken  at  a  height  0.1  in.  above  the  base  of 
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Fig.  5-1  Arrangement  for  Spectroscopic  Measurements 
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the  plasma  jet.  During  a  traverse  the  voltage  drop  across  tie  plasma 
torch  was  recorded  every  fifteen  seconds.  The  net  power  input  to  the 
plasma  was  calculated  for  each  voltage  drop  reading  and  these  values 
were  averaged  to  compute  the  average  net  power  for  each  traverse. 

The  radiation  seen  by  the  spectroscope  within  any  bandwidth  is 
the  sum  of  the  spectral  line  radiation  and.  the  continuum  radiation 
within  that  bandwidth.  In  order  to  determine  the  intensity  of  the 
spectral  line,  two  traverses  were  made  at  each  power  level.  The  first 
traverse  was  made  with  the  spectroscope  set  to  receive  all  of  the  radia¬ 
tion  underneath  the  spectral  line,  including' the  continuum.  The  data 
thus  obtained  is  referred  to  as  the  total  intensity.  The  spectroscope 
was  then  set  just  far  enough  away  from  this  frequency  range  to  exclude 
the  spectral  line  radiation,  thus  measuring  just  the  continuum  intensity. 
The  continuum  intensity  was  subtracted  from  the  total  intensity  to  yield 
the  lateral  spectral  line  intensity.  The  details  of  the  settings  of 
the  spectroscope  for  these  various  measurements  are  discussed  in  Appendix 
E. 

The  two  pairs  of  traces  which  were  used  to  determine  the  tempera¬ 
ture  distribution  in  the  plasma  jet  are  shown  in  Figures  5-2,  5-5*  5-4* 
and  5~5*  The  chart  speed  was  4  in/min  so  that  one  inch  on  the  recorder 
chart  corresponds  to  0.025  in.  in  the  plasma  jet.  The  argon  flow  rate 
to  the  plasma  torch  in  these  traverses  was  5*1  lbs/hr;  the  cooling 
water  flow  rate  was  1,260  Ibs/hr.  The  average  net  power  input  for 
the  traverses  in  Figures  5-2  and  5-5  was  10,100  Btu/hr,  for  Figures 
5-4  and  5~5>  15*000  Btu/hr. 

The  data  were  taken  from  a  smooth  curve  drawn  directly  .on  the 
spectroscope  recorder  traces.  The  maximum  of  each  curve  was  determined 
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14.5  14.0  13.5  13.0  12.5  12.0  11.5 

TIME  (min.) 


Fig.  5-3  Traverse  of  Plasma  Jet;  Continuum 
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Fig.  5-4  Traverse  of  Plasma  Jet;  Argon  Atom  Line  4158.59A  and  Continuum 
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Fig.  5-5  Traverse  of  Plasma  Jet;  Continuum 
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and  a  vertical  line  was  dropped  from  the  maximum.  The  relative  intensity 
was  read  at  horizontal  increments  of  0.2  in.  from  the  maximum  on  the 
recorder  trace,  which  corresponds  to  horizontal  increments  of  0.005  in* 
in  the  plasma  jet.  The  values  of  the  relative  intensity  at  equal  dis¬ 
tances  on  both  sides  of  the  maximum  were  averaged.  This  averaging  pro¬ 
cedure  amounts  to  folding  the  recorder  trace  along  a  vertical  line  at 
the  maximum  and  drawing  an  average  curve  through  the  two  superposed 
curves.  The  folding  process  was  used  to  adjust  for  any  asymmetry  in 
the  plasma  jet  and  to  eliminate  the  asymmetry  introduced  by  the  hysteresis 
in  the  recorder.  This  procedure  was  used  on  both  the  total  and  continuum 
intensity  traces.  The  folded  continuum  intensity  was  then  subtracted 
from  the  folded  total  intensity  at  corresponding  horizontal  increments 
to  yield  the  relative  lateral  spectral  line  intensity  at  the  two  power 
levels.  These  data  are  shown  in  Figure  5-6  and  5-7. 

Nestor  and  Olsen  (Ref.  18)  have  developed  a  numerical  technique 
based  on  the  Abel  integral  equation  to  invert  the  lateral  intensity 
distribution  to  a  radial  intensity  distribution.  The  occurrence  of  the 
Abel  integral  equation  in  this  connection  is  discussed  in  Appendix  E. 

This  numerical  technique  has  been  programmed  for  use  on  an  I.B.M.  709 
digital  computer  by  Enright  (Ref.  17).  The  lateral  spectral  line. in¬ 
tensities  shorn  in  Figures  5-6  and  5-7  were  inverted  to  radial  spectral 
line  intensities  using  the  computer  program.  The  computer  output  was 
normalized  at  its  maximum  and  plotted.  The  plots  of  the  radial  spectral 
line  intensity  are  shown  in  Figures  5-8  and  5-9®  The  maxima  on  these 
curves  corresponds  to  a  temperature  of  l5?550°K  in  the  plasma  jet.  The 
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RELATIVE  INTENSITY 
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° 

Fig.  5-9  Radial  Intensity;  Argon  Atom  Line  4158.59A;  Average  Net  Power  15,000  Btu/hr 
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temperature  at  horizontal  increments  of  O.OOlj.  in.  from  the  axis  of  the 
plasma  jet  was  determined  by  comparing  the  radial  spectral  line  intensity 
with  the  theoretical  relative  intensity  of  the  spectral  line  (see 
Appendix  E).  The  temperature  profiles  in  the  argon  plasma  jet  deter¬ 
mined  at  the  two  power  levels  are  shown  in  Figure  £-10. 

5.3  DETERMINATION  OF  RADIAL  DISTRIBUTION  OF  RADIATIVE  POWER  DENSITY  IN  AN 
ARGON  PLASMA  JET 

The  bolometer  and  plasma  torch  were  positioned  as  shown  in  Figure 
£-11  with  the  optical  system  providing  a  five-fold  magnification  of  the 
image  of  the  plasma  jet  in  the  plane  of  the  bolometer  sensing  element. 

The  bolometer  was  calibrated  in  this  arrangement  as  discussed  in  Appen¬ 
dix  B.  Care  was  taken  to  place  the  bolometer  element  at  a  position  in 
the  image  of  the  plasma  jet  corresponding  to  the  height  in  the  plasma  jet 
at  which  the  lateral  spectral  line  intensity  was  measured. 

The  lateral  intensity  distribution  of  the  continuum  radiation  was 
measured  by  traversing  the  bolometer  horizontally  across  the  image  of 
the  plasma  jet.  Readings  of  the  bolometer  output  signal  were  recorded 
at  0.0£0  in.  increments  at  the  lower  intensities,  and  at  0„02£  in. 
increments  at  the  higher  intensities.  These  increments  correspond  to 
0.010  in.  and  0.00£  in.  increments  in  the  plasma  jet,  respectively. 
Readings  of  the  voltage  drop  across  the  plasma  torch  were  recorded 
simultaneously  with  each  bolometer  output  signal.  The  net  power  input 
to  the  plasma  was  calculated  for  each  reading.  These  values  were  then 
averaged  to  compute  the  average  net  power  for  each  traverse. 

The  lateral  intensity  distribution  of  the  continuum  radiation 
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Fig.  5-10  Temperature  Profiles  of  Plasma  Jet 
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was  measured  at  three  power  levels®  The  data  from  these  measurements 
are  shown  in  Figures  $-12,  $-13,  and  5-lU.  The  smooth  curves  drawn 
through  these  points  were  folded  at  the  maximum  intensity  and  a  symmetri¬ 
cal  lateral  intensity  distribution  was  drawn  between  the  experimental 
curves  for  each  traverse.  These  folded  curves  are  shown  in  Figures  $- 1 5, 
5-16,  and  5-17.  The  lateral  intensity  of  the  continuum  radiation  was 
read  from  these  folded  curves  at  O.OOli  in.  increments.  These  data  were 
inverted  to  radial  distributions  using  the  computer  program  (Ref.  17). 

The  resulting  radial  distribution  of  the  radiative  power  density  in 
the  argon  plasma  jet  at  three  power  levels  is  shown  in  Figure  5-l8. 

5.4  DETERMINATION  OF  CONTINUUM  RADIATION  AS  A  FUNCTION  OF  TEMPERATURE 

Ideally,  the  lateral  intensity  of  the  continuum  radiation  should 
have  been  measured  at  the  same  power  levels  as  those  at  which  the 
temperature  was  determined.  However,  it  is  nearly  impossible  to  operate 
the  plasma  torch  at  a  predetermined  power  level,  because  the  voltage 
drop  across  the  torch  cannot  be  precisely  predicted.  As  a  result,,  the 
radiation  traverses  were  taken  at  different  power  levels  than  the 
temperature  traverses.  A  linear  interpolation  technique  was  used  to 
adjust  the  temperature  profiles  and  radial  distributions  of  the  radia¬ 
tive  power  density  to  the  same  power  levels.  The  method  of  interpola¬ 
tion  and  sample  calculations  are  shown  in  Appendix  F.  Figure  $-.19  shows 
the  final  results  of  the  temperature-radiation  determination.  The 
theoretical  models  discussed  in  Ghapter  k  are  shown  for  comparison. 
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Fig.  5-12  Lateral  Radiative  Flux;  Average  Net  Power  12,100  Btu/hr 
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RADIATIVE  FLUX  (watls/cm2) 


Fig.  5-14  Lateral  Radiative  Flux;  Average  Net  Power  17,200  Btu/hr 
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Fig.  5-15  Folded  Lateral  Radiative  Flax 
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Fig.  5-16  Folded  Lateral  Radiative  Flux 
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Fig.  5-17  Folded  Lateral  Radiative  Flux 
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Fig.  5-18  Radiant  Energy  Transfer  as  a  Function  of  Radius 
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Fig.  5-19  Radiation  from  an  Atmospheric  Argon  Plasma  as  a  Function  of  Temperature 
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5.5  DISCUSSION  OF  RESULTS 

About  three  minutes  is  required  to  make  one  traverse  across  the 
image  of  the  plasma  jet  with  the  spectroscope.  Since  two  traverses  are 
required  to  determine  the  lateral  spectral  line  intensity  distribution, 
the  plasma  torch  must  operate  with  negligible  power  fluctuations  for 
about  seven  minutes.  Due  to  the  fluctuations  of  the  voltage  drop 
across  the  plasma  torch,  it  was  very  difficult  to  maintain . steady  opera¬ 
tion  of  the  torch  for  this  length  of  time.  In  practice,  as  many  as  ten 
pairs  of  traverses  were  made  before  one  pair  was  obtained  during  which 
the  power  remained  nearly  constant.  Notice  on  Figure  5-2  the  sudden 
jumps  in  the  recorder  chart  trace  between  8,25  min,  and  8,50  min.  This 
is  an  example  of  the  difficulties  caused  by  variations  in  the  voltage 
drop  across  the  plasma  torch.  In  the  case  of  this  particular  trace,  the 
voltage  was  unsteady  for  only  about  ten  seconds  and  then  returned  to  its 
original  value,  A  smooth  curve  was  drawn  through  the  unsteady  region 
and  the  data  were  taken  from  this  curve.  If  the  voltage  had  settled  to 
a  new  value  following  this  disturbance,  the  traverse  would  have  been 
useless. 

The  spectroscopic  traces,  Figures  5-2,  5-3,  5-U,  and  5-5,  show  an 
asymmetry  about  the  maximum.  Some  of  this  asymmetry  is  due  to  hysteresis 
in  the  recorder,  but  some  is  also  caused  by  an  asymmetry  in  the  plasma 
jet.  The  latter  effect  is  due  to  some  misalignment  of  the  nozzle  over 
the  electrode,  which  causes  the  arc  to  be  distributed  non-uniforraly 
around  the  electrode.  The  effects  of  the  asymmetry  do  not  appear  in 
the  final  results  because  of  the  folding  procedure  which  was  used.  The 
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alignment  of  the  nozzle  was  found  to  be  very  critical  in  reducing  this 
asymmetry. 

The  data  of  the  lateral  radiative  flux.  Figures  5-12,  5-13*  and 
5-lil,  show  some  scatter  in  the  data  points  and  an  asymmetry  about  the 
maximum.  The  scatter  is  due  to  voltage  fluctuations.  As  about  thirty 
minutes  was  required  to  make  each  of  these  traverses,  it  was  impossible 
to  avoid  this  scatter.  The  error  introduced  because  of  this  scatter 
and  the  asymmetry  shows  up  in  Figure  5-18.  Although  it  is  difficult  to 
determine  the  magnitude  of  this  error,  it  does  not  appear  to  be  prohibi¬ 
tive,  judging  from  the  results. 

The  final  results  shown  in  Figure  5-19  also  exhibit  scatter  in 
the  data  points.  This  scatter  is  not  severe  except  in  the  points  above 
19,000°K,  This  scatter  and  any  error  in  the  results  are  the  accumulation 
of  experimental  errors  and  the  errors  incurred  in  the  manipulation  of 
the  data.  The  experimental  curve  has  been  drawn  by  hand  as  the  best 
apparent  representation  of  the  results.  Due  to  the  scatter  of  the 
points  above  19,000°K,  this  portion  of  the  curve  cannot  be  considered 
as  reliable  as  the  portion  below  19,000°K. 

The  experimental  results  indicate  that  both  bremmstrahlung  and 
recombination  are  important  sources  of  the  radiation  emitted  from  an 
atmospheric  argon  plasma  over  the  range  of  temperatures  investigated. 

The  Kramers -Uns Bid  model  predicts  the  radiation  emitted  from  an  atmos¬ 
pheric  argon  plasma  reasonably  well  from  10,000°K  to  17,000°K. 

There  are  several  possible  reasons  for  the  discrepancy  between 
the  experimental  data  and  the  Kramers-UnsOld  model.  The  extensive 
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experimental  measurements  and.  data  processing  necessary  to  arrive  at 
the  results  increases  the  possibility  of  errors.  Some  of  these  errors 
show  up  as  scatter  in  the  data,  hut  a  systematic  error  could  contri¬ 
bute  to  the  discrepancy  between  the  theoretical  model  and  the  experi¬ 
mental  data.  The  cutoff  frequency,  which  is  not  accurately  known, 
has  been  assumed  to  be  independent  of  temperature.  Finally,  the 
approach  used  of  adding  the  radiation  predicted  by  the  Kramers  model 
from  the  doubly  ionized  atoms  to  that  predicted  from  the  singly  ionized 
atoms  is  probably  not  a  good  approximation. 

A  model,  shown  in  Figure  5-20,  was  formulated  to  check  the  fold¬ 
ing  technique  and  method  of  inversion  used  to  obtain  i  he  experimental 
curve  on  Figure  5-19*  The  section  through  the  thickest  portion  of  the 
plasma  jet  was  divided  into  rectangular  elements  of  volume  0.01  cm. 
thick  having  an  area  of  one  square  centimeter.  The  value  of  the  radia¬ 
tion  at  each  radius  was  taken  from  the  experimental  curve  on  Figure 
5-19*  The  contributions  from  the  elements  were  added  and  the  result¬ 
ing  value  was  compared  with  the  maximum  lateral  radiative  flux  mea¬ 
sured.  The  calculations  are  shown  in  Appendix  F.  The  agreement  with 
the  measured  values  is  good,  which  verifies  the  general  method,  of 
arriving  at  the  results. 
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CHAPTER  6 
CONCLUSIONS 

The  first  conclusion  to  he  drawn  from  the  experimental  results  is 
that  radiation  is  an  important  mode  of  heat  transfer  from  an  argon 
plasma  jet  generated  in  a  plasma  torch  and  exhausted  to  atmospheric 
pressure.  Radiant  energy  transfer  from  the  argon  plasma  jet  studied 
accounts  for  as  much  as  18  per  cent  of  the  net  power  input  to  the  plasma. 

The  determination  of  the  actual  radiative  heat  transfer  from  the 
argon  plasma  and  the  plasma  torch  indicates  that  a  large  amount  of 
energy  is  radiated  from  the  core  of  plasma  within  the  exit  nozzle. 

The  intensity  of  the  radiation  from  the  hot  core  of  plasma  greatly 
exceeds  the  intensity  of  the  radiation  from  the  plasma  jet.  This  addi¬ 
tional  radiation  is  significant  and  must  be  considered  in  an  overall 
energy  balance  about  a  plasma  torch. 

The  energy  radiated  from  an  atmospheric  nitrogen  plasma  jet  is 
less  than  1  percent  of  the  net  power  input  to  the  plasma,  up  to  100,000 
Btu/hr  net  power  input.  The  total  radiation  from  a  nitrogen  plasma  jet 
exhausted  to  the  atmosphere  is  negligible  for  the  flow  rate  and  nozzle 
conf iguration  investigated. 

The  radiation  from  an  atmospheric  argon  plasma  which  is  assumed 
to  be  in  local  thermodynamic  equilibrium  can  be  predicted  reasonably 
well  with  the  Kramers-L'nsflld  model  from  10,000°K  to  17,000°K. 

The  thermistor  bolometer  is  well  suited  to  measure  the  radiation 
emitted  by  a  laboratory  plasma.  It  is  sensitive  over  a  wide  range  of 
the  spectrum  and  its  calibration  curve  is  linear.  The  measurements 
are  made  without  disturbing  the  flow  in  the  plasma  jet. 
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APPENDIX  A 

EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

A-1  PLASMA  GENERATING  FACILITY 

A  schematic  layout  of  the  plasma  generating  facility  in  the  North¬ 
western  Gas  Dynamics  Laboratory  is  shown  in  Figure  A-1.  Direct  current 
power  to  operate  the  plasma  torch  is  supplied  by  four  Model  SRH  I4I4I4. 

Miller  Electric  Company  arc  welders.  Each  welder  is  capable  of  producing 
about  ten  kilowatts  of  d.c.  power.  The  welders  are  arranged  so  that 
ary  number  of  units  can  be  used  to  supply  power  to  the  torch.  The  units 
may  be  connected  in  series,  in  parallel,  or  in  a  series-parallel  combina-. 
tion.  A  filter  circuit  is  available  which  reduces  the  a.c.  ripple  in 
the  welder  output  to  less  than  one  per  cent  of  the  nominal  current. 

No  change  in  the  radiation  from  the  plasma  was  observed  in  any  of  the 
experimental  measurements  upon  incorporating  the  filters  into  the  power 
circuit.  The  leads  from  the  welders  to  the  torch  carry  both  the  current 
and  the  cooling  water.  A  high  frequency  starter  manufactured  by  Thermal 
Dynamics  Corporation  is  used  to  provide  the  breakdown  voltage  to 
initiate  the  arc.  The  welders  are  controlled  by  varying  the  voltage  on 
the  control  windings  of  the  transformers  in  the  welders. 

Water  for  cooling  the  torch  is  taken  from  the  building  supply. 

An  auxiliary  pump  is  available  to  increase  the  flow  rate  of  the  cooling 
water  above  that  obtained  using  just  tap  pressure.  This  pump  was  not 
used  in  this  experimental  work.  Nominal  control  of  the  water  flow  rate 
is  provided  by  partially  closing  the  shut-off  valve  from  the  building 
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supply. 

The  gas  used  to  generate  the  plasma  is  supplied  from  ordinary 
high  pressure  cylinders.  Gases  of  the  highest  purity  available  commer¬ 
cially  were  purchased  from  a  local  distributor  of  the  Linde  Company. 

The  entire  plasma  generating  facility  is  operated  from  a  remote 
control  panel.  Controls  for  the  power  supply — welders,  filter  circuit, 
and  high  frequency  starter — as  well  as  various  gauges,  flow  meters,  and 
recorders  are  housed  in  the  control  panel. 

The  d.c.  current  to  the  torch  was  determined  by  measuring  the 
potential  across  a  1^00  ampere  $0  millivolt  shunt  manufactured  by  the 
Weston  Electric  Instrument  Company  which  is  in  series  with  the  power 
supply.  This  potential  was  measured  with  a  Varian  Associates,  Inc., 

Model  G  11  A  strip  chart  recorder.  The  input  resistance  of  this  recor¬ 
der  has  a  minimum  value  of  l/2  megohm.  The  recorder  was  calibrated 
before  each  series  of  experimental  runs  using  a  Leeds  and  Northrup  No. 
8662  portable  potentiometer.  A  typical  recorder  trace  made  during  one 
of  the  experimental  runs  is  shown  in  Figure  A-2. 

The  voltage  drop  across  the  torch  was  measured  with  a  Model  14.30 
Weston  voltmeter.  Clip  leads  from  this  voltmeter  were  connected  directly 
across  the  torch.  With  this  arrangement,  the  voltmeter  measures  only 
the  voltage  drop  across  the  plasma  torch  and  does  not  include  the 
voltage  drop  in  the  electrical  leads  from  the  welders  to  the  torch. 

Determination  of  the  power  loss  from  the  torch  to  the  cooling 
water  was  accomplished  by  measuring  the  flow  rate  of  the  cooling  water 
and  its  temperature  rise.  The  flow  rate  was  measured  with  a  Fischer 
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Porter  No,  B5-27-5/77  rotameter  located  in  the  control  panel.  The 
temperatures  of  the  incoming  and  outgoing  cooling  water  were  measured 
with  copper-cons tantan  thermocouples.  The  thermocouples,  mounted  in 
capillary  glass  tubing  sealed  with  wax,  were  inserted  in  the  water 
pipes  directly  upstream  and  downstream  of  the  torch.  The  constantan 
leads  from  the  thermocouples  were  commoned  so  that  the  e.m.f.  between 
the  copper  leads  was  proportional  to  the  temperature  rise  of  the  cool¬ 
ing  water.  This  e.m.f.  was  measured  with  a  Leeds  and  Northrup  Type  G 
Speedomax  recording  potentiometer.  A  typical  trace  made  by  the  thermo¬ 
couple  recorder  during  one  of  the  experimental  runs  is  shown  in  Figure 
A-3. 

The  flow  rate  of  the  gas  used  to  generate  the  plasma  was  measured 
with  a  Fischer  Porter  No.  2F-1/H-25-5/70  Tri-Flat  variable  area  flow 
meter  located  in  the  control  panel.  A  stainless  steel  float  was  used 
in  the  flow  meter.  The  pressure  and  temperature  of  the  incoming  gas 
were  indicated  by  gauges  mounted  in  the  control  panel. 

The  seeder  used  to  seed  the  plasma  jet  is  a  Sylco  CCC  Mark  VH 
fine  powder  feeder  manufactured  by  Sylvester  and  Company,  In  the  seeder, 
the  powder  is  metered  from  a  vibrating  storage  tank  by  a  feed  screw  to 
a  carburetor  where  it  is  dispersed  in  the  carrier  gas  stream.  The 
feed  screw  has  a  variable  speed  drive  which  allows  continuous  adjustment 
of  the  powder  feed  rate  from  zero  up  to  the  maximum  feed  rate.  A  #57x51 
nozzle  insert  manufactured  by  Thermal  Dynamics  Corporation  was  used  in 
the  plasma  torch  when  seeding  the  plasma  jet.  This  nozzle  has  a  small 
hole  through  which  the  carrier  gas  and  powder  are  injected  radially 
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Fig.  A-3  Typical  Temperature  Recorder  Plot 
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into  the  plasma  jet. 

The  carbon  particles  with  which  the  plasma  jet  was  seeded  are 
manufactured  by  Cabot  Corporation  under  the  name  of  Supercarbovar. 

O 

The  particles  have  a  mean  diameter  of  l.lj.  x  10  m.  (Ref.  I;).  The  same 
high  purity  argon  used  to  generate  the  plasma  was  used  as  a  carrier  gas 
through  the  seeder. 

A-2  EXPERIMENTAL  PROCEDURE 

The  following  start-up  procedure  was  used  at  the  beginning  of  each 
experimental  run.  The  bolometer  power  supplies  were  turned  on  and, 
after  a  warm-up  period  of  about  one  minute,  the  noise  level  in  the  bolo¬ 
meter  circuit  was  read  on  the  output  voltmeter.  No  experimental  measure¬ 
ments  were  made  when  the  noise  level  exceeded  three  millivolts.  The 
electrical  current  recorder  and  the  thermocouple  recorder  were  positioned 
with  the  recording  pens  on  a  unit  interval  chart  line,  and  the  charts 
were  titled  with  a  run  number.  The  current  recorder  was  calibrated  at 
the  beginning  of  each  series  of  experimental  runs.  Since  the  thermo¬ 
couple  recorder  has  an  internal  calibration  circuit  it  required  only 
occasional  calibration. 

The  gas  supply  to  the  torch  was  turned  on  and  the  pressure  was 
adjusted  with  the  gas  regulator  to  80  p.s.i.g.  The  cooling  water  was 
turned  on  and  the  flow  rate  was  adjusted  to  the  desired  value.  A  flow 
rate  was  chosen  which  would  give  a  maximum  temperature  rise  of  about 
25°F. 

The  gas  flow  rate  through  the  torch  was  set  at  its  maximum  value 
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of  about  12  lbs/hr „  With  the  main  power  to  the  torch  turned  offs  the 
high  frequency  starter  was  held  on.  The  ring  of  spark  produced  by  the 
starter  in  the  annular  region  between  the  electrode  and  nozzle  was  ob¬ 
served  through  the  exit  nozzle  in  a  mirror  held  above  the  torch.  The 
position  of  the  nozzle  over  the  electrode  was  adjusted  so  that  the  spark 
was  evenly  distributed  around  the  electrode.  The  gas  stream  issuing 
from  the  exit  nozzle  was  checked  visually  for  water  droplets  which  would 
indicate  a  water  leak  in  the  torch.  In  the  absence  of  a  water  leak* 
the  main  power  to  the  torch  was  turned  on.  The  arc  was  struck  by 
pressing  the  starter  switch  momentarily.  The  gas  flow  rate  was  then 
reduced  to  the  desired  operating  value. 

With  the  torch  in  operation,,  the  clip  leads  from  the  Weston  volt¬ 
meter  were  connected  across  the  torch®  The  current  recorder  and  a  stop¬ 
watch  were  started  simultaneously.  One  minute  later  the  thermocouple 
recorder  was  started.  Alter  a  final  check  of  all  of  the  instruments 3 
the  experimental  readings  were  begun. 

The  gas  temperature 9  pres sure 3  and  flow  rate ,  and  the  cooling 
water  flow  rate  were  recorded  at  the  beginning  and  end  of  each  experi¬ 
mental  run.  The  gas  pressure  and  flow  rate  were  observed  to  remain 
constant  during  the  runs.  The  gas  temperature  was  observed  to  drop 
very  gradually  during  the  runs.  This  is  due  to  the  cooling  accompanying 
the  expansion  of  the  gas  in  the  high  pressure  cylinder.  This  tempera¬ 
ture  drop  was  never  more  than  3°F  during  an  experimental  run.  The 
cooling  water  flow  rate  was  observed  to  fluctuate  randomly  +  3  per  cent. 
The  values  of  the  gas  pressure  and  flow  rate  used  in  the  calculations 
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were  those  recorded  at  the  beginning  of  each  run.  The  average  of  the 
gas  temperature  at  the  beginning  and  end  of  each  run  was  used  in  the 
calculations.  The  mean  value  about  which  the  flow  rate  of  the  cooling 
water  fluctuated  was  used  in  the  calculations.  This  mean  value  did  not 
vary  during  an  experimental  run. 
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APPENDIX  B 

CALIBRATION  OF  BOLOMETER 

B-l  CALIBRATION  APPARATUS 

The  bolometer  was  calibrated  using  a  black  surface  of  known  tempera¬ 
ture.  The  calibration  for  the  total  radiation  measurements  was  accom¬ 
plished  by  heating  a  hollow  copper  cone  in  a  Hevi-Duti  Type  62  electric 
furnace.  The  cone,  shown  in  Figure  B-l,  was  located  so  that  its  interior 
surface  was  visible  from  outside  the  furnace  through  a  hole  in  the  end 
of  the  furnace.  Due  to  repeated  heating  and  cooling,  a  heayy  scale  of 
copper  oxide  was  formed  on  the  entire  surface  of  the  cone,  rendering  it 
a  nearly  black  surface.  Iron-constantan  thermocouples  were  located 
along  the  length  of  the  cone  in  such  a  way  that  the  temperature  of  the 
area  of  the  cone  viewed  by  the  bolometer  was  known.  A  radiation  shield 
made  of  asbestos  board  served  to  define  the  area  from  which  the  radia¬ 
tion  was  emitted  and  to  shield  the  bolometer  from  the  stray  radiation 
emitted  by  the  warm  exterior  surface  of  the  furnace.  A  vertical  copper 
plate  was  used  as  a  radiation  source  to  calibrate  the  bolometer  for  use 
with  the  optical  systems  which  are  described  below.  Its  temperature  was 
also  measured  with  an  iron-constantan  thermocouple.  The  vertical  copper 
plate  served  to  define  a  focal  plane  for  the  optical  systems.  The 
thermocouple  potentials  were  measured  with  the  Leeds  and  Northrup  No. 

8662  portable  potentiometer.  The  electric  furnace  has  a  maximum 
operating  temperature  of  1850°F, 
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Fig.  B-l  Experimental  Set-Up  for  Total 


Radiation  Calibration  (vertical  scale 


magnified  8  fold) 
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B-2  CALIBRATION  OF  BOLOMETER  FOR  TOTAL  RADIATION  MEASUREMENTS 

Consider  an  element  of  black  surface ,  dA^,  radiating  in  all 
directions  from  one  side.  Some  of  this  radiation  is  intercepted  by 
black  element,  dAg,  at  a  distance  r  from  dA^.  The  line  connecting  these 
two  surface  elements  forms  angles  <£  ^  and  with  the  normals  to  dA-^ 
and  dA^,  respectively.  The  rate  of  radiant  energy  transfer,  dq^  g,  from 
dA-^  to  dA^  is  (Ref.  19) 


i-2 


Ii  o'Ai  cos  <f>\  dAz  cos  <f>z 


(B-1-) 


1^  is  the  intensity  of  radiation  from  element  dA^.  The  elemental  solid 
angle  subtended  by  dAg  at  dA^  is 


du\  ~ 


dfKz  COS  4>z 


Substituting,  Eq,  B-l  becomes 

dq{_z  =  X,  dk |  cos  <j>{  (b-2) 

The  rate  of  radiant  energy  transfer,  dq^,  from  dA^  in  all  directions  is 
found  by  integrating  over  the  complete  hemispherical  angle  above  dA-^ 

35T  *  x\  /cos  h  d“\  (B"3) 

i  0 

The  quantity  on  the  left  hand  side  of  Eq.  B-3  is,  by  definition,  the 
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emissive  power  of  the  black  surface  dA^ 
4  _  d(\\ 


W,  -  o- V  = 


o' A . 


where  <r  is  the  Stefan-Boltzmann  constant,  0.1713  x  10“®  Btu/hr  ft^  °R^ 
is  the  temperature  of  black  surface  A-^  (°R) 

The  quantity  on  the  right  hand  side  of  Eq.  B-3  may  be  integrated  by 
noting  that,  due  to  the  symmetry  of  the  problem. 


dWy  3= 


Zirr  sin  </>,  r  dfy 


3=  2 t  sin  d<f> | 


Substituting  these  quantities  into  Eq.  B-3 


tt/Z 


<r  T, 


X,  J  2r  cos  sin  d<f>y 


L  ■  2, 


,7T/2 


=  Ij  y  sin2<£ 


=  Ii7r 

Thus  the  proportionality  constant  in  Eq.  B-l  for  a  black  surface  element, 
dA^,  is 


i,  =  fV 


Eq.  B-l  may  now  be  used  to  evaluate  the  radiant  heat  transferred 
through  an  imaginary  surface,  Ag,  at  a  distance  r  from  a  black  surface 
A^.  For  the  purpose  of  calibrating  the  bolometer,  the  black  surface, 
A^,  was  taken  as  the  hole  in  the  radiation  shield  in  front  of  the  hot 
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black  surface  of  the  interior  of  the  cone.  The  object  is  to  find  the 
radiative  flux,  F,  at  the  bolometer  element  at  distance  r  from  the  hole 
in  the  radiation  shield.  Since  the  bolometer  was  positioned  on  the  axis 
of  the  cone,  and  since  the  plane  of  the  radiation  shield  was  perpendicu¬ 
lar  to  the  axis  of  the  cone, 

=  0®,  <f>2  £  0® 

For  a  distance  r  from  the  hole  to  the  bolometer  element  much  greater 
than  the  radius,  r^,  of  the  hole  in  the  radiation  shield 


*,  =  0® 

Likewise,  since  the  bolometer  element  is  very  small  (0.5  mm  x  0.5  mm) 

*2  =  0® 

Eq,  B-l  becomes 


d<k  i-2 


o-T,4  tfA,  d A2 
v  r2 


(B-h) 


bolometer  element  due  to  the 
hole,  Eq,  B-l;  is  integrated 


However,  since  the  bolometer  element  is  very  small,  dAg>  may  be  taken 
as  the  area,  SAg,  of  the  bolometer  element.  Performing  the  integration 


In  order  to  find  the  radiative  flux  at  the 
radiation  from  the  entire  area,  A1 ,  of  the 
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Bizi  =  <r1i4Ai 

8A2  t rr2 


(B-5) 


The  quantity  on  the  left  of  Eq.  B-f>  is  recognized  as  the  radiative  flux, 
F,  at  the  bolometer  element.  Thus,  the  calibration  equation  is 


F  = 


a  T|  A, 


irr 


(B-6) 


2 

where  F  is  the  radiative  flux  at  the  bolometer  element  (Btu/hr  ft' ) 

o-  is  the  Stefan-Boltzmann  constant,  0.1713  x  10"®  Btu/hr  ft^  °R^ 
T^  is  the  temperature  of  the  black  surface  (°R) 
is  the  area  of  the  black  surface  (ft  ) 
r  is  the  distanoe  from  the  bolometer  element  to  the  black  surface 


(ft) 

As  an  example,  one  of  the  calibration  points  taken  supplied  the 
following  data: 

temperature  of  surface  of  cone  in  view  of  bolometer,  1,23S>°F 
distance  from  bolometer  element  to  radiation  shield,  2k  in. 
diameter  of  hole  in  radiation  shield,  21/32  in. 
bolometer  output  signal,  19*8  mv 

Using  Eq.  B-6  to  find  the  flux  at  the  bolometer  element, 


(0.1713  X  I0~8)(I235  +  460)V)(l/4)(0.656/12)2 
rr(  24/12  )2 

*  2.64  Btu/hr  ft2 
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To  find  the  calibration  value,  the  flux  at  the  bolometer  element  is 
divided  by  the  corresponding  bolometer  output  signal  in  volts,  giving 

=  133  /hr  ft*  volt 

Figure  B-2  shows  the  calibration  data  which  were  taken  by  Dowdy 
and  Bowman."  The  data  were  taken  with  the  bolometer  at  three  positions 
along  the  axis  of  the  cone.  The  temperature  of  the  cone  was  detected 
by  three  thermocouples,  located  as  shown  in  Figure  B-l.  The  scatter 
of  the  points  on  Figure  B-2  is  due  to  the  fact  that  the  temperature  of 
the  cone  was  not  constant  along  its  length.  Referring  to  Figure  B-l, 
the  area  of  the  cone  viewed  by  the  bolometer  in  position  1  is  predomi¬ 
nantly  in  the  vicinity  of  thermocouple  1,  while  the  area  of  the  cone 
viewed  by  the  bolometer  in  positions  2  and  3  is  predominantly  in  the 
vicinity  of  thermocouple  2.  Using  only  these  data  points,  a  straight 
line  was  drawn  through  the  points  as  shown  in  Figure  B-3.  This  is  the 
calibration  curve  which  was  used  for  all  of  the  total  radiation  measure¬ 
ments.  This  value  for  the  calibration  was  checked  several  times  over 
the  months  during  which  the  experimental  data  were  taken.  The  calibra¬ 
tion  value  was  always  found  to  be  within  one  per  cent  of  that  originally 
determined  as  shown  in  Figure  B-3.  The  sample  calculation  of  a  cali¬ 
bration  value  shown  above  was  taken  from  one  of  the  later  calibration 
checks. 

*>£ 

"Research  Assistants,  Department  of  Mechanical  Engineering,  Northwestern 
University. 
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Fig.  B-3  Calibration  Curve  for  Total  Radiation 
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B-3  CALIBRATION  OF  BOLOMETER  FOR  USE  WITH  OPTICAL  SYSTEMS 

The  calibration  of  the  bolometer  for  use  with  the  optical  systems 
shown  in  Figures  3->7  and  5-H  was  somewhat  simpler  than  the  calibration 
technique  used  for  total  radiation  measurements*  For  calibration  with 
the  optical  systems ,  the  copper  plate  was  used  as  a  radiation  source  in 
the  furnace®  The  furnace  and  copper  plate  were  positioned  such  that  the 
plane  of  the  copper  plate  defined  one  focal  plane  of  the  optical  system® 
With  the  bolometer  element  in  the  other  focal  plane  of  the  optical 
system,  the  bolometer  element  is  effectively  placed  directly  on  the 
surface  of  the  copper  plate®  Due  to  this  fact,  it  is  unnecessary  to 
take  into  account  the  angular  distribution  of  the  radiation  from  the 
black  surface  as  was  done  in  Section  B.2.  Considering  the  copper  plate 
as  a  black  surface  at  temperature  T,  the  radiative  flux  at  the  surface 
of  this  copper  plate  is  o-T^4.  Dividing  the  value  of  the  radiative  flux 
at  the  surface  of  the  plate  by  the  corresponding  output  signal  of  the 
bolometer  gives  the  calibration  value.  This  calibration  value  already 
includes  any  effects  of  absorption  within  the  optical  system. 

As  an  example,  the  data  from  one  of  the  calibration  points  using 
the  5:1  magnification  optical  system  shown  in  Figure  5-11  were  as  follows? 
temperature  of  copper  plate,' 1,588°F 
bolometer  output  signal,  31  mv 

The  radiative  flux  at  the  surface  of  the  copper  plate  is 

crT4  s  (O.I7I3  x  I0”®)(I588  4-  460)4  =  2.96  x  I04  Bfu/hr  ft2 

The  calibration  value  is  .found  by  dividing  the  radiative  flux  at  the 
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surface  of  the  plate  by  the  bolometer  output  signal  in  volts 

2.96  X  10*  Btu/hr  ft2  =  Q  96  x  (0©  Btu/hr  ft2  volt 
3.1  X  ICT2  volts 

The  calibration  curves  obtained  for  use  with  the  optical  systems 
are  shown  in  Figures  B-ii  and  B-5. 
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RADIATIVE  FLUX  AT  SOURCE  (Btu/hr  ft*x  I0'4) 


Fig.  B-4  Calibration  of  Bolometer  for  Pointwise  Measurements  (2:1  magnification) 
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BOLOMETER  OUTPUT  SIGNAL  (volts) 


Fig.  B-5  Calibration  of  Bolometer  for  Pointwise  Measurements  (5:1  magnification) 
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APPENDIX  C 

SAMPLE  CALCULATIONS  OF  RESULTS  IN  CHAPTER  3 

C-l  CALCULATION  OF  TOTAL  RADIATION  FROM  PLASMA  JET 

One  of  the  total  radiation  points  shown  on  Figure  3-3  was  calcu¬ 
lated  from  the  following  data: 
distance  from  bolometer  to  plasma  jet  I4O  3/8  in. 
voltage  drop  across  plasma  torch  21.8  v 
cooling  water  flow  rate  823  lbs/hr 
current  to  torch  3^8  amps 
cooling  water  temperature  rise  l3«7°F 
bolometer  output  signal  77  mv 
The  gross  power  input  to  the  torch  is 

(348  amps)(2l.8  volts)  =  7,590  watts 

(7,590  watts )(3.4I5  Btu/watt  hr)  =  25,900  Btu/hr 

The  power  loss  to  the  cooling  water  is 

(825  lbs/hr)(I.O  Btu/lb  °F)(I5.7  °F)  =  12,900  Btu/hr 

The  net  power  input  to  the  plasma  is 


25,900  Btu/hr  -  12,900  Btu/hr  =  13,000  Btu/hr 

The  area  of  the  sphere  on  which  the  bolometer  element  lies  is 

4 it  (  40.6  m  \2  _  ^2 

M2  in/ft'  TT 
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The  radiative  flux  at  the  bolometer  element  is 

(0.077  volts)(  133  Btu/hr  ft2  volt)  *  10.2  Btu/hr  ft2 

The  total  power  radiated  from  the  plasma  jet  (radiant  power  crossing 
the  sphere  in  Figure  3-1)  is 


(144  ft2 )(I0.2  Btu/hr  ft2)  =  1,470  Btu/hr 


The  per  cant  net  power  radiated  is 


1,470  Btu/hr 
13,000  Btu/hr 


(100%) 


11.3% 


C-2  CALCULATION  OF  TOTAL  RADIATION  FROM  POINTWISE  DATA 

One  of  the  experimental  runs  from  which  pointwise  radiation  data 
was  taken  was  made  under  the  following  conditions: 
optical  system,  two-fold  magnification 
average  net  power  input,  11,300  Btu/hr 
height  of  image  of  plasma  jet,  2  l/8  in® 

The  diameter  of  the  base  of  the  cone,  Figure  3-9 ,  was  taken  as  3/l6  in., 
the  diameter  of  the  exit  nozzle  opening.  Radiation  readings  were  taken 
at  1/8  in.  vertical  increments  along  the  axis  of  the  image  of  the  plasma 
jet.  As  an  example,  consider  the  reading  which  was  taken  at  a  height 
of  1/2  in.  along  the  centerline  of  the  plasma  jet.  The  area  of  a 
frustum  of  a  right  circular  cone  having  a  radius  r^  of  the  base  and 
of  the  top  and  an  altitude  h  is 
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,r<ri  +  r2>  X  V/hS+ir|-r2>2 

Considering  a  vertical  section  through  the  axis  of  the  conical  plasma 
jet,  the  data  for  this  example  are  found  by  considering  the  similar 
triangles 

5/32  =  r, 

17/16  17/16  -  |5/32 

r8  =  0.087  in 

5/32  s  r2 

17/16  17/16  -  17/32 

r2  =  0.078  in 

The  height  of  the  frustum  under  consideration  in  the  plasma  jet  is 
l/l6  in.  The  area  of  this  frustum  is 

tt(.087  +  .078)v/(.063)2  +  (.087  -  078)2  =  0.033  in2 

The  bolometer  output  signal  at  this  point  was  3.9  volts  (see  Figure  3-8). 
Using  the  calibration  value  from  Figure  B-U,  the  radiative  flux  at  the 
surface  of  the  plasma  jet  at  the  vertical  midpoint  of  this  frustum  of 
the  cone  is 

(3.9  voits)(l.94  x  I05  Btu/hr  ft2  volt)  =  7.6  X  I05  Btu/hr  ft2 

The  radiant  power  crossing  the  surface  of  this  frustum  of  the  cone  is 
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,0.033  in2  v 

(  144  jn2/ft2  M7  6  x  l0*  Btu/hr  ft2)  *  173  Btu/hr 

The  total  radiation  from  the  cone  was  calculated  by  adding  the 
contribution  from  each  frustum  of  the  cone,  calculated  as  in  the  example 
above*  The  average  net  power  input  to  the  plasma  was  calculated  by 
averaging  the  net  power  input  at  each  radiation  reading.  The  per  cent 
net  power  radiated  from  the  plasma  jet  was  calculated  by  dividing  the 
total  radiation  from  the  cone  by  the  average  net  power  input  to  the 
plasma. 

C-3  CALCULATION  OF  ACTUAL  RADIATIVE  HEAT  TRANSFER 

The  actual  radiative  heat  transfer  from  the  plasma  and  torch  was 
calculated  using  the  data  on  the  angular  distribution  of  the  radiation 
from  the  torch  by  imagining  the  bolometer  element  to  lie  on  a  sphere 
whose  center  is  the  center  of  rotation  of  the  plasma  jet.  The  data  were 
taken  with  the  bolometer  lt7  5/8  in.  from  the  center  of  rotation.  As  an 
example,  consider  the  radiation  reading  taken  at  a  viewing  angle  of  20°. 
The  readings  were  taken  at  5°  increments  at  the  higher  intensities. 

The  reading  at  20°  is  considered  to  lie  in  the  center  of  the  area  of  the 
sphere  included  between  the  viewing  angles  of  17.5°  and  22, 5°.  The 
area  of  a  spherical  segment  of  height  h  of  a  sphere  of  radius  r  is 

2ir  rh 

In  terms  of  the  quantities  known  in  the  experiment 
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h  *  r(l  -  cos  9) 

where  9  is  the  viewing  angle#  The  area  of  the  spherical  segment 
between  the  viewing  angles  of  0°  and  17*5°  is 

g^fn/fV8*1  ~ cos  l7S°>  =  455  f,z 

The  area  of  the  segment  of  the  sphere  between  the  viewing  angles  of  0° 
and  22o5°  is 

Zv  (^27fn/ft)2d  "  cos  22.5°)  =  7.52  ft2 

The  area  of  the  ring-like  element  of  the  sphere,  across  which  the  radia¬ 
tion  measured  at  a  viewing  angle  of  20°  passes,  is 

7.52  ft2  -  4.55  ft2  =  2.97  ft2 

Consider  the  radiation  reading  taken  at  a  viewing  angle  of  20° 
with  the  torch  operating  at  an  average  net  power  input  of  8,700  Btu/hr. 
The  operating  conditions  at  this  point  were : 
voltage  drop  across  torch  23#0  v 
cooling  water  flow  rate  6kB  lbs/hr 
current  to  torch  2lj0  amps 
cooling  water  temperature  rise  15>*f>°F 
bolometer  output  signal  110  rnv 
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The  gross  power  input  to  the  torch  is 

(23.0  volts)(240  amps)(3.4l5  Btu/watt  hr)  =  18,800  Btu/hr 

The  power  loss  to  the  cooling  water  is 

(645  lbs/hr)(I.O  Btu/lb  °F)(I5.5  °F)  =  10,000  Btu/hr 

The  net  power  input  to  the  torch  is 

18,800  Btu/hr  -  10,000  Btu/hr  =  8,800  Btu/hr 

The  radiative  flux  at  this  viewing  is 

(0.110  volts) (133  Btu/hr  ft2  volt)  *  14.6  Btu/hr  ft2 

The  contribution  to  the  actual  radiative  heat  transfer  passing  across 
the  segment  of  the  sphere  between  17.5°  and  22.5°  viewing  angle  is 

(2.97  ft2)(l4.6  Btu/hr  ft2)  =  43.4  Btu/hr 

The  actual  radiative  heat  transfer  from  the  plasma  and  torch  was 
calculated  by  adding  the  contributions  passing  across  each  ring-like 
segment  of  the  sphere.  The  average  net  power  input  to  the  plasma  was 
calculated  by  averaging  the  net  power  input  at  each  radiation  reading. 
The  per  cent  net  power  radiated  from  the  plasma  and  torch  was  calculated 
by  dividing  the  total  radiative  power  by  the  average  net  power  input. 
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APPENDIX  D 

SAMPLE  CALCULATIONS  FOR  THEORETICAL  MODELS 

D-l  SAMPLE  CALCULATION  FOR  KRAMERS  APPROXIMATION 

The  Kramers  approximation  is  given  by 

P(T)  =  1.4  X I0~34  Z2!^  T,/2  [^»-84°/T  -  *-53,800/1 j  watts/em3  ()i-7 ) 

As  an  example,  consider  a  temperature  of  19,000°K.  From  Drellishak's 
tables  (Ref,  10)  for  the  eo^ilibrium  composition  of  atmospheric  argon 
plasma 

n^  =  1.81;  x  1017  cm-3 
n'F  =  1.0  8  x  lO"*"3  cm""3 

l 

n  =  1.87  x  1017  cm-3 

e 

I  II 

where  n. ,  n.  are  the  number  densities  of  the  first  and  second  ion, 
x  * 

respectively.  The  contribution  to  the  radiation  from  free-free  electron- 
first  ion  collisions  is 

Pj-CT)  =  (1.4  X  I0“34)(|)2(l.87  X  I0i7)(I.84XI0i7)(I.9X|04),/£  x 

f _ ! _ _  _  _i_  1 

X  |_*.0442  ff2.83  J 

PX(T)  =  596  watts/cm3 

The  contribution  to  the  radiation  from  free-free  electron-second  ion 
collisions  is 
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£(T)  =  (1.4  X  I0-34)(2)2(I.87  X  I0l7)(l .08  X  I0I8)(I.9X  I04)1'*  X 

X  [  £.0442  "  ff2.83  ]  =  14  Watts/cm3 

The  radiation  predicted  by  the  Kramers  approximation  at  19,000°K  is 

o  o  o 

3 96  watts/cnr  +  II4  watts/cnr  =  610  watts/ crrr 

D-2  SAMPLE  CALCULATION  FOR  SOMMERFELD  APPROXIMATION 

The  Sommerfeld  approximation  is  given  by 


P(T)  =  l.4X|0"34Z2nenjTl/2g(T)[tf"840/T-tf“63'80°/T]  watts/cm3  (ii-8) 

Again,  consider  a  temperature  of  19,000°K.  The  average  Gaunt  factor, 
g(T),  at  19,000°K  is  given  by  Greene  (Ref.  11)  as  1.31.  The  radiation 
predicted  at  19,000°K  by  the  Sommerfeld  approximation  is 
(610  watts/cm3)  (1.31)  =  800  watts/cm3 

D-3  SAMPLE  CALCULATIONS  FOR  KRAMERS-UNSOLD  APPROXIMATION 

The  emission  coefficient  given  by  the  Unsdld  approximation  is 

<sv  =  6.8  x  I0“46  Z2  watt  sec /cm3  (k-9) 

The  emission  coefficient  given  by  the  Kramers  approximation  i.s 

=  6.8  X  I0“45  Z2  e-'nv/WT  WQtf  sec/Cm3  &-3) 
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The  combined  Kramers-UnsiJld  model  forces  the  emission  coefficient  of 
the  Kramers  approximation  to  match  the  emission  coefficient  of  the 
Unsold  approximation  at  the  cutoff  frequency,  vq  ,  by  adjusting  the 
constant  term  in  the  Kramers  emission  coefficient*  Equating  the  two 
emission  coefficients  at  the  cutoff  frequency, 

68  X  I0“45  Z2  =  CZ2  neni  e-hvq/kT 

jl/2  jl/2 

which  gives 

C  =  6.8  x  |0"45  e  hvg/kT  (d-i) 

as  the  constant  term  in  the  Kramers  emission  coefficient.  The  contribu¬ 
tion  to  the  radiation  below  the  cutoff  frequency  from  free-free  and 
free-bound  transitions  is  found  by  integrating  the  Unsold  emission  co¬ 
efficient  from  the  lower  limit  of  transmissibility  of  the  sodium  chloride 
window,  P|  ,  to  the  cutoff  frequency,  . 

P(T)  =  /  9€V  dv  =  6.8  X  icr45  Z2^(v  -*,)  (D-2) 

v\  T  * 

The  contribution  to  the  radiation  above  the  cutoff  frequency  from  free- 
free  transitions  is  found  by  integrating  the  Kramers  emission  coeffi¬ 
cient  from  the  cutoff  frequency,  ,  to  the  upper  limit  of  trans¬ 
missibility  of  the  sodium  chloride  window,  vu 


113 


AEDC-TDR-64-233 


P<T) 


/“cz2  fil  dv 

4  T1'2 

CZ2  (k/h)  nen|  T1'2  (e-hxgAT  -  h*u/!<T  ) 


Substituting  for  C  and  evaluating  the  constant  term  gives 


P(T)  =  1.42  x  I0“34  e  l4.400/T  z2  n@nj  j,/2  x 


X^-14,400/7  -  ^-53,800/j) 


(D-3) 


The  upper  and  lower  frequencies  of  transmissibility  of  the  sodium 
chloride  window  are 
vu  «  11*2  x  10114  sec"1 
vj  =  0.175  x  10 ^  sec"1 

The  cutoff  frequency  is  taken  as 
Vq  -  3.0  x  10  sec 

and  is  assumed  to  be  constant  with  temperature. 

Consider,  as  an  example,  a  temperature  of  19,000°K.  The  contri¬ 
bution  to  the  radiation  predicted  by  the  Unsold  approximation  is 


P(T) 


(6.8  X  I0"46)(l)2 


(1.87  X  I0i7)(I,84  X  IQ17) 
(1.9  x  I04)1'2 


(2.82  X  |014) 


478  watts /cm3 
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The  contribution  to  the  radiation  predicted  by  the  modified  Kramers 
approximation  is 


P(T)  =  (1.42  x  I0"34)(ff758)(l.87x  IO,7)[(l)2(l.84  x|017)  + 

+  (2)z(l.08  x  I015)]  (1.9  X  l04)l/2(*--758  -  *~283) 
=  602  watts  /cm3 

The  radiation  predicted  by  the  Kramers-UnsBld  model  at  19,000°K  is 
1*78  watts/cm^  +  602  watts/cm^  =  1,080  watts/cm^ 
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APPENDIX  E 

TEMPERATURE  DETERMINATION 

E-l  THEORY  OF  TEMPERATURE  DETERMINATION 

The  electrons  populating  the  higher  energy  levels  of  the  excited 
atoms  and  ions  in  a  plasma  emit  radiation  of  characteristic  frequencies 
upon  undergoing  transitions  to  lower  energy  levels.  The  intensity  of 
the  radiation  of  any  spectral  line  is  given  by  (Ref.  16) 

_  -  a  h  n(T)ga  *"Ea /kT 

Xafi  ~  Aa£  ”*q/9  (E-l) 

2  9j  tf-Ei/kT 

i  =  l 

o 

where  Iayg  is  the  intensity  of  the  radiation  (watts/crrr) 

Aa0  is  the  transition  probability  from  state  e  to  state 
frecluency  Of  the  radiation  (sec-'1') 
n(T)  is  the  number  density  of  the  atoms  or  ions  emitting  the 
radiation  (cm  ■3) 

"t  h 

g^  is  the  statistical  weight  of  the  i  state 

k  is  the  Boltzmann's  constant,  1.38  x  10“  ergs/  K 

T  is  the  temperature  (°K) 

E^  is  the  excitation  energy  of  the  ith  state 

The  frequency  of  the  spectral  line  is  given  by  Planck's  equation 

t>»«0  =  Ea-  Ep  (E-2) 

The  quantity  in  the  numerator  of  Eq.  E-l  is  the  number  atoms 
(considering  an  atom  spectral  line)  in  the  a  state  as  given  by  the 
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Boltzmann  distribution.  The  quantity  in  the  denominator  of  Eq.  E-l  is 
the  internal  partition  function  of  the  atom.  As  the  temperature  rises 
at  constant  pressure ,  the  decrease  in  the  density  of  the  gas  and  the 
destruction  of  atoms  by  ionization  overtakes  the  increasing  exponential 
and  the  intensity  exhibits  a  maximum. 

The  theoretical  relative  intensity  of  the  i|l58.59  5  argon  atom 
line.  Figure  E-l,  was  known  (Ref#  13).  The  maximum  intensity  of  this 
line  occurs  at  l5,330°K.  if  the  temperature  at  some  point  in  the  plasma 
jet  is  above  15>,££0°K,  the  experimentally  measured  spectral  line  intensity 
will  also  show  a  maximum#  In  this  case  it  is  unnecessary  to  measure  the 
absolute  intensity  of  the  spectral  line;  the  temperature  in  the  plasma 
jet  can  be  determined  from  the  relative  intensity  of  the  spectral  line 
by  comparison  with  the  theoretical  relative  intensity  of  the  line. 

The  intensity  of  radiation  emitted  by  a  non-opaque  plasma  depends 
upon  the  depth  of  the  plasma  emitting  the  radiation.  For  the  purpose 
of  determining  the  temperature  distribution  in  the  plasma  jet,  the 
plasma  is  assumed  to  be  optically  thin.  This  assumption  is  supported 
by  measurements  in  an  argon  plasma  jet  (Ref.  13).  Under  this  assump¬ 
tion,  the  intensity  of  the  characteristic  radiation  of  a  spectral  line 
emitted  by  the  plasma  depends  on  the  temperature  distribution  in  the 
plasma  jet  and  the  thickness  of  the  jet  at  the  point  of  observation. 

The  plasma  jet  is  also  assumed  to  have  axial  symmetry.  Thus,  the 
cross-section  at  any  vertical  height  in  the  plasma  jet  is  circular. 

The  temperature  distribution  at  any  height  in  an  axially  symmetrical 
plasma  jet  will  be  a  function  only  of  the  radius. 
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Fig.  E-l  Theoretical  Relative  Intensity;  Argon  Atom  Line  4158.59A 
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E-2  ABEL  INTEGRAL  EQUATION 

The  intensity  of  the  radiation  viewed  by  the  spectroscope  is  the 
sura  of  radiation  emitted  by  each  volume  element  throughout  the  depth 
of  the  plasma.  Figures  5-6  and  5-7  show  the  spectral  line  intensity  as 
functions  of  the  lateral  distance  from  the  axis  of  symmetry  of  the  plasma 
jet.  The  lateral  spectral  line  intensity  is  an  intensity ,  integrated 
over  the  depth  of  the  plasma. 

Assuming  axial  symmetry  within  the  plasma  jet,  each  horizontal 
cross-section  is  a  circle.  Define  a  rectangular  coordinate  system  on 
this  circle  with  its  origin  at  the  center  of  the  circle  as  shown  in 
Figure  E-2  such  that  x  represents  the  lateral  distance  from  the  axis  of 
syrmnetry  at  which  the  spectroscope  observes  the  radiation  from  the  plasma. 
Denote  the  lateral  spectral  line  intensity  shown  in  Figures  5-6  and  5-7 
by  I(x).  The  object  is  to  find  a  relationship  between  I(x)  and  l(r), 
the  spectral  line  intensity  as  a  function  of  the  radial  coordinate  with¬ 
in  the  circle.  The  lateral  spectral  line  intensity  is  given  in  terms 
of  the  radial  spectral  line  intensity  by  the  integral 


I(x) 


I(r)  dy 


By  definition  of  the  coordinate  systems 


r2  =  x2  +  y2 

Y  is  the  value  of  the  coordinate  y  on  the  boundary  of  the  circular 
cross-section  of  radius  it,  i.e., 
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Fig.  E-2  Derivation  of  Abel  Integral 
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R2  -  x2  +  Y2 

Changing  variables  in  the  integral  at  a  constant  x 


xu)  -  2  xiiist 

This  is  the  Abel  integral  equation,  which  can  be  solved  analytically 
(Ref.  20)  giving 


X(r) 


I  ,R  X'U)  d% 

wI 


where  I'(x)  is  the  derivative  of  I(x)  with  respect  to  x 
Olsen  (Ref.  16)  has  developed  a  numerical  technique  to  evaluate  this 
integral  using  data  of  the  form  of  that  in  Figures  $-6  and  £-7.  Enright 
(Ref.  17)  has  programmed  this  technique  for  use  on  an  I.B.M.  709  digi¬ 
tal  computer.  Enright's  program  was  used  to  invert  the  experimental 
data. 


E-3  SPECTROMETER  ADJUSTMENTS  FOR  THE  MEASUREMENT  OF  LATERAL 
SPECTRAL  LINE  INTENSITY 

The  spectrum  of  the  radiation  from  an  argon  plasma  consists  of  a 
series  of  spikes  superposed  on  the  slowly  varying  continuum  radiation. 
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The  spikes  correspond  to  various  spectral  lines  of  the  atoms  and  ions. 
Figure  E-3  shows  a  spectrum  of  the  plasma  taken  in  the  course  of  these 
experiments.  This  spectrum  was  taken  with  a  minimum  entrance  and  exit 
slit  opening  to  provide  maximum  resolution. 

Theoretically,  the  spectral  lines  should  be  very  narrow.  However, 
the  spectral  lines  in  a  plasma  are  broadened  due  to  various  temperature 
and  electric  field  effects.  In  order  to  measure  the  intensity  of  any 
spectral  line,  it  is  necessary  to  open  the  exit  slit  of  the  spectrometer 
so  as  to  include  all  of  the  radiation  under  the  desired  spectral  line. 

Ideally,  a  spectral  line  should  be  chosen  for  the  temperature 
determination  which  is  isolated  from  all  other  lines.  However,  this  is 
not  the  only  factor  to  be  considered  in  the  choice.  The  Ul£8 .59  X  atom 
line  was  chosen  by  other  investigators  for  several  reasons  (Refs.  13,  16, 
17).  The  hl$Q &9  X  line  is  flanked  by  two  other  lines,  one  at  Ul5U»50  X 
and  one  at  I4I6I4..I8  X.  With  the  exit  slit  opened  to  receive  all  of  the 
radiation  under  the  U158.59  X  line,  the  radiation  from  these  other  two 
lines  will  be  included,  Enright  (Ref.  17)  estimates  that  the  U15U.50  X 
line  contributes  a  negligible  amount  of  radiation  to  the  readings.  He 
also  states  that,  although  the  radiation  from  the  iO.6U.l8  X  line  is  not 
negligible,  the  variation  of  the  intensity  of  this  line  with  temperature 
is  nearly  the  same  as  for  the  10-58 .59  X.  Therefore,  a  negligible  error 
is  introduced  by  including  these  two  lines  in  the  spectroscopic  measure¬ 
ment  of  the  lateral  intensity  of  the  U158.59  X  line  for  the  purpose  of 
determining  the  temperature  distribution  in  the  plasma  jet. 

The  necessary  exit  slit  width  was  determined  by  measuring  the  width 
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of  the  U158.59  S,  line  on  Figure  E-3  and  dividing  this  width  by  the 
linear  dispersion  of  the  spectrometer  at  this  wavelength,  h3  S/ran 
(Ref®  17)®  The  exit  slit  width  used  was  208  microns  greater  than  the 
entrance  slit  width®  A  spectrum  of  the  plasma  taken  with  the  exit  slit 
widened  is  shown  in  Figure  E— 14 .  The  loss  of  resolution  due  to  the 
wider  exit  slit  opening  is  pronounced. 

From  the  spectrum  in  Figure  E-I4,  the  setting  of  the  spectrometer 
corresponding  to  the  maximum  intensity  in  the  region  of  the  ljl58.5>9  $. 
line  was  determined.  The  measurements  of  the  total  lateral  intensity 
of  the  line  and  continuum  were  made  at  this  spectrometer  setting.  In 
order  to  determine  what  portion  of  this  total  lateral  intensity  was  due 
to  continuum  radiation,  the  spectrometer  was  adjusted  to  receive  the 
continuum  radiation  just  to  the  left  of  the  ltl58 ,£9  &  line  in  Figure 
E-U,  but  excluding  all  of  the  radiation  from  the  line®  A  second  traverse 
was  made  at  this  spectrometer  setting  to  measure  the  continuum  intensity. 
The  continuum  intensity  was  assumed  to  be  the  same  underneath  the  itlS>8.59 
X  line  as  slightly  to  the  left  of  it,  i.e®,  at  slightly  shorter  wave¬ 
lengths,  The  continuum  intensity  was  measured  at  a  wavelength  of  1*135  X. 

The  traverses  were  made  with  an  entrance  slit  width  of  10  microns 
and  an  exit  slit  width  of  218  microns.  The  entrance  slit  was  opened 
this  additional  amount  to  increase  the  intensity  of  radiation  received 
by  the  photomultiplier  tube*  The  resulting  additional  loss  of  resolu¬ 
tion  does  not  effect  the  temperature  determination. 
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appendix  f 

CALCULATIONS  OF  TEMPERATURE-RADIATION  DATA 

F-l  INTERPOLATION  OF  TEMPERATURE  RADIATION  DATA  TO  COMMON  POWER  LEVELS 
Ideally,  the  temperature  profiles  of  the  plasma  jet,  Figure  5-10, 
and  the  inverted  radiation  profiles.  Figure  5-18,  should  have  been 
determined  at  the  same  average  net  power  input  to  the  plasma.  However, 
as  has  been  explained  in  Section  5.U,  it  is  impossible  to  operate  the 
plasma  torch  at  a  predetermined  power  level*  As  a  result,  two  tempera¬ 
ture  profiles  were  determined  at  an  average  net  power  input  of  10,100 
Btu/hr  and  15,000  Btu/hr,  while  the  three  radiation  profiles  were  deter' 
mined  at  an  average  net  power  input  of  12,100  Btu/hr,  13,900  Btu/hr 
and  17,200  Btu/hr.  In  order  to  determine  the  radiation  as  a  function 
of  temperature,  these  data  were  adjusted  by  means  of  a  linear  inter¬ 
polation  technique  to  common  power  levels. 

To  calculate  the  points  on  Figure  5-19  corresponding  to  12,100 
Btu/hr  average  net  power,  the  two  temperature  profiles  were  used  to 
interpolate  to  a  temperature  profile  at  12,100  Btu/hr.  At  each  O.OOlt 
in®  along  the  radius,  the  temperature  difference,  AT(r),  between  the 
temperature  profile  at  15,000  Btu/hr  and  the  temperature  profile  at 
10,100  Btu/hr  was  calculated  from  the  profiles  on  Figure  5-10*  The 
temperature  difference  AT’(*0  between  a  hypothetical  temperature  pro¬ 
file  at  12,100  Btu/hr  and  the  measured  temperature  profile  at  10,100 
Btu/hr  was  calculated  at  each  0o00U  in®  along  the  radius  using  the 
interpolation  technique  of  the  following  formula 
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/ 12,100  -  10,100 

'  15,000  -  10,100 


AT'(r) 


The  temperature  at  each  O.OOl;  in.  along  the  radius  corresponding  to  an 
average  net  power  input  of  12,100  Btu/hr  was  found  by  adding  the  tempera¬ 
ture  difference,  AT’(r),  to  the  temperature  at  the  same  radius  measured 
at  10,100  Btu/hr,  following  the  formula 

Tl 2,100  *  TIO,IOO*r>  +  AT'(r> 

Using  this  method,  the  temperature  profile  measured  at  10,100  Btu/hr 
can  be  interpolated  to  any  desired  power  level,  P^,  between  10,100 
Btu/hr  and  13,000  Btu/hr.  The  interpolation  formula  is 


Tp  (r) 

1) 


TI0,I00  (r> 


+  (  %  ~  10,100  \  AT  (r) 

'  15,000  -  10,100/ 


(F-l) 


This  equation  was  used  to  interpolate  to  temperature  profiles  at  12,100 
Btu/hr  and  13,900  Btu/hr  average  net  power  input.  With  the  temperature 
profile  known  at  the  same  power  levels  as  each  radiation  profile,  the 
points  on  Figure  3-19  were  plotted. 

In  order  to  avoid  extrapolating  the  temperature  profile  past  the 
measured  values,  a  similar  interpolation  technique  was  used  to  find  a 
radiation  profile  at  13,000  Btu/hr  by  using  the  radiation  profiles  at 
13,900  Btu/hr  and  17,200  Btu/hr.  The  interpolation  formula  used  is 


Rl5,000(r) 


Rl3,900*r) 


+  )  AR(r)  (F-2) 

N  17,200  -  13,900' 
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where  AR(r)  is  the  difference  in  radiation  at  each  radius  r  between 
the  radiation  profile  at  17,200  Btu/hr  and  the  radiation 
profile  at  13 , 900  Btu/hr „ 

The  radiation  profile  at  15,000  Btu/hr  obtained  by  applying  Eq.  F-2  was 
compared  with  the  measured  temperature  profile  at  15,000  Btu/hr  to  deter¬ 
mine  the  third  set  of  points  on  Figure  5-19. 

For  example,  consider  a  radius  of  0.020  in.  in  the  plasma  jet. 

From  Figure  5-10,  the  temperature  at  a  radius  of  0.020  in.  at  15,000 
Btu/hr  is  22,500OKj  the  temperature  at  10,100  Btu/hr  is  18,050°K. 

Thus  the  temperature  difference,  AT(r),  at  a  radius  of  0.O2O  in.  is 
3,950°K.  Using  Eq.  F-l,  the  temperature  corresponding  to  an  average 
net  power  input  of  12,100°K  at  this  radius  is 


TIO  ,««(  020)  =  18,050  +  ( 12,100  ~  l0,l0°  )( 3,950) 

12,100  x  15,000”  10,100  ' 

=  19,660  ®K 

Referring  to  Figure  5-18,  the  radiation  at  a  radius  of  0.020  in.  at  an 

O 

average  net  power  input  of  12,100  Btu/hr  is  1,380  watts/cm  «  This  is 
one  of  the  points  plotted  on  Figure  5=19 •  Likewise,  at  a  radius  of 
0.020  in®  the  temperature  corresponding  to  an  average  net  power  input  of 
13,900  Btu/hr  is 


/ 13,900  -  10,100  \ 

T|3>9Oo  (-020)  -  18,050  +  fe^oO^JOO  ) {3,950) 
=  21,110  °K 
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Referring  to  Figure  5-18,  the  radiation  at  a  radius  of  0,020  in,  at  an 
average  net  power  input  of  13,900  Btu/hr  is  1,200  watts/cm  ,  This  is 
another  point  on  Figure  5-19® 

Again,  at  a  radius  of  0,020  in,,  the  radiation  from  Figure  5-18 

•a 

at  an  average  net  power,  input  of  17,200  Btu/hr  is  l,it70  watts/cm  ,  and 
the  radiation  at  13,900  Btu/hr  is  1,220  watts/cm^.  The  difference  in 
radiation  from  these  two  profiles  at  a  radius  of  0,020  in,  is  260 
watts/cm  ,  Using  Eq»  F-2,  the  radiation  corresponding  to  an  average 
net  power  input  of  15,000  Btu/hr  at  this  radius  is 


RI5.000<  020'  -  '•22°  +  (Ir’zoo-  13,900  )(260) 


=  1,310  watts  /cm3 

From  Figure  5-10,  the  temperature  at  a  radius  of  0.020  in,  and  at  an 
average  net  power  input  of  15,000  Btu/hr  is  18,300°K,  This  is  also  a 
point  on  Figure  5-19. 

F-2  CHECK  ON  TEMPERATURE-RADIATION  DATA 

A  check  on  the  experimental  curve  drawn  on  Figure  5-19  was  made. 
Using  the  interpolation  technique  discussed  above,  the  temperature  at 
0.00U  in.  increments  along  the  radius  was  found  at  an  average  net  power 
input  of  12,100  Btu/hr.  The  plasma  was  imagined  to  be  divided  into 
incremental  volumes  one  centimeter  square  and  0,01  centimeters  thick 
along  a  diameter  as  shown  in  Figure  5-20.  The  radiation,  in  watts/cm  , 
was  found  at  each  increment  along  the  radius  for  an  average  net  power 
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input  of  12,100  Btu/hr,  by  using  the  temperature  at  each  increment  as 
calculated  with  the  interpolation  formula,  Eq„  F-l,  and  reading  the 
radiation  corresponding  to  each  temperature  from  the  experimental  curve 
on  Figure  5-19«  The  radiation  emitted  by  each  elemental  volume  in 
Figure  5-20  is  0.01  times  the  value  read  from  Figure  5-19?  since  the 
volume  of  each  element  in  Figure  5-20  is  0.01  cm  *  The  contributions 
to  the  radiation  from  the  elemental  volumes  along  the  radius  were  added 
and  the  result  was  doubled  to  include  the  contributions  along  the 
diameter  in  the  plasma  jet.  Doubling  the  contributions  from  the  ele¬ 
ments  along  the  radius  counts  the  elemental  volume  at  the  center  of  the 
plasma  jet  twice.  Therefore,  the  contribution  from  the  elemental  volume 
at  the  center  was  subtracted  once.  The  resulting  figure  should  agree 
with  the  maximum  value  of  the  lateral  radiative  flux  which  was  measured 
at  12,100  Btu/hr,  Perforating  the  calculation  indicated  gives  a  maximum 
lateral  radiative  flux  of  kOh  watts/cm2  .  Referring  to  Figure  5-12, 
the  maximum  measured  value  of  the  radiative  flux  was  1+18  watts/cm2  .  A 
similar  calculation  of  the  radiative  flux  at  an  average  net  power  input 
of  13,900  Btu/hr  gives  a  value  of  1*28  watts/cm2  .  Referring  to  Figure 
5-13,  the  maximum  measured  value  of  the  radiative  flux  at  this  power 
level  was  1*26  watts/cm2  . 

Referring  to  Figure  5-lU,  the  maximum  measured  radiative  flux  at 
an  average  net  power  input  of  17,200  Btu/hr  was  lj.65  watts/ cm2  .  This 
value  may  be  interpolated  to  correspond  to  an  average  net  power  input 
of  15,000  Btu/hr  in  a  manner  similar  to  that  indicated  in  Eq.  F-2. 

The  calculation  gives 
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max 


/  15,000  -  13,900  \ 

(15,000)  =  426  +  (  — -) (465-  426) 

1.0, ww/  v  17,200  -  13,900  7 


=  439  watts /cm1 


Using  the  temperature  profile  at  an  average  net  power  input  of  15,000 
Btu/hr  on  Figure  5-10,  and  the  experimental  curve  on  Figure  5-19,  the 
calculated  maximum  radiative  flux  is  lj.50  watts/cnr  . 

These  results  may  be  summarized  as  follows. 


Average 

Net  Power 

Maximum  Lateral  Radiative  Flux 

%  Error 

Input 

watts/ 

cur 

(from 

measured 

(Btu/hr) 

Calculated 

Measured 

value) 

12,100 

ltOU 

1118 

-2.9 

13,900 

U28 

126 

+o.5 

15,000 

U5o 

U39 

+2.5 
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